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Abstract 

We present a formalism for the calculation of multi-particle one-loop amplitudes, valid 
for an arbitrary number of external legs, and for massive as well as massless particles. 
A new method for the tensor reduction is suggested which naturally isolates infrared di- 
vergences by construction. We prove that for N > 5, higher dimensional integrals can 
be avoided. We derive many useful relations which allow for algebraic simplifications of 
one-loop amplitudes. We introduce a form factor representation of tensor integrals which 
contains no inverse Gram determinants by choosing a convenient set of basis integrals. 
For the evaluation of these basis integrals we propose two methods: An evaluation based 
on the analytical representation, which is fast and accurate away from exceptional kine- 
matical configurations, and a robust numerical one, based on multi-dimensional contour 
deformation. The formalism can be implemented straightforwardly into a computer pro- 
gram to calculate next-to-leading order corrections to multi-particle processes in a largely 
automated way. 
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1 Introduction 



The quest for new physics at TeV coUiders hke the Fermilab Tevatron, the upcoming Large 
Hadron Colhder at CERN, as well as the International Linear Collider (ILC) project, requires 
the quantitative calculation of many hard multi-particle processes. Direct searches rely on the 
proper identification of multi-particle signals and a precise understanding and determination 
of the corresponding multi-particle/multi-jet backgrounds. The latter, in particular those for 
multi-jet or vector boson(s) plus multi-jet production, are considerable, as can be estimated 
from leading order studies ^ |2] . Indirect searches at high luminosity machines will also involve 
high precision measurements of multi-particle observables of the Standard Model, comprising 
six-point processes like e^e~ —>■ 4 fermions at the ILC, so that accurate predictions for these 
cross sections will be mandatory. High energy physics is thus entering a new era where the 
quantitative description of hard multi-particle production is becoming a topic of prominent 
importance, whereas the lowest order estimates for such processes are plagued by the well- 
known deficiencies of large renormalisation and factorisation scale dependences, poor multi-jet 
modelling and large sensibility to kinematic cuts. Therefore the calculation of next-to-leading- 
order (NLO) corrections to such cross sections is a necessary step forward. However, this task 
involves an enormous technical complexity. 

To perform an NLO calculation with hard external particles, the following points have 
to be addressed: 

1. Generation of tree amplitudes with + 1 external particles 

2. Subtraction of soft and collinear singularities due to single unresolved real radiation 

3. Generation of the one- loop amplitude with external particles 

4. Evaluation of the loop diagrams, UV renormalisation, extraction of soft and collinear 
singularities 

5. Combination of the contributions above, cancellation of soft singularities, cancellation of 
collinear singularities or absorption into distribution functions 

6. Numerical evaluation of the finite amplitude 

For step 1, efficient, highly automated tools and algorithms are available. A similar state- 
ment can be made for step 2. Although automatisation is less trivial for this point, the avail- 
able algorithms [HI IH El E] are well tested and in principle valid for an arbitrary number of 
legs. The same is not true for loop amplitudes. Although efficient programs like FeynArts [7j, 
GRACE ISl or QGRAF ilOj exist which reliably deal with the combinatorial complexity of 
generating multi-leg one-loop Feynman diagrams, the evaluation of these diagrams, in partic- 
ular in the presence of infrared divergences and for more than four external legs, is still far 
from being automated. So the bottleneck in constructing an automated program package for 
NLO cross sections is step 3, the evaluation of the loop graphs. Although a number of five- 
parton processes, see for example have been calculated already, and very recently even 

physical 2 — >^ 4 results have become available ISS] . these calculations all required a tedious 
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individual treatment, and most of them do not allow to obtain fully differential results. There- 
fore it is desirable to have tools which allow the calculation of NLO cross sections in a largely 
automated way. These tools should be able to handle massless as well as massive particles, and 
should be numerically reliable and fast. Ideally, they should also allow to be interfaced with a 
parton shower in a universal way, using for example the formalism proposed in or jSSl ■ 

Several approaches to this aim have been suggested in the literature so far, from purely 
numerical ones to ones where the emphasis is on algebraic manipulations. A completely nu- 
merical approach has been worked out by D. Soper et al. jSElEZl; where the sum over cuts for a 
given graph is performed before the numerical integration over the loop momenta. In this way 
unitarity is exploited to cancel soft and coUinear divergences before they show up as explicit 
poles. 

However, the conventional method of calculating the virtual (loop) and real (radiation) parts 
separately, thus generating infrared poles which cancel in the sum, is still the most widely and 
successfully used approach so far. Of course, within this approach, there are still many different 
ways to proceed, in particular in what concerns the evaluation of the one-loop amplitude. The 
most straightforward procedure - and historically the first one - relies on the use of recursion 
relations to reduce the tensor integrals occurring in the one-loop amplitude to a set of known 
basis integrals jSH]-|HZI- In the recent work on this subject, the emphasis is primarily on methods 
which are suitable for an efficient numerical evaluation of multi-leg amplitudes. For the massless 
case, a formalism has been proposed recently in j33], which produces spurious inverse Gram 
determinants, but in [S3] a method is proposed how to deal with them. The formalism given 
in j3T] avoids inverse Gram determinants in the reduction of pentagon integrals, but deals with 
massive particles only. In j^ni; another algorithm is presented, using spinor helicity methods. 
Based on the formalism of an evaluation of one- loop integrals in massless gauge theories 
for up to 12 external legs has been given recently in [57^. 

A numerical approach to the one-loop integrals is the one of [HE], where various concepts 
like the Bernstein-Tkachov theorem, Mellin-Barnes representation and sector decomposition 
are combined to get a stable numerical behaviour in all regions of configuration space. A fully 
numerical approach to the calculation of loop integrals by contour integration also has been 
elaborated in [HH] ■ A semi-numerical approach, where a subtraction formalism for the UV and 
soft /coUinear divergences of the one-loop graphs has been worked out, is presented in jUU]. The 
idea is to integrate the remaining finite part in loop momentum space without performing any 
tensor reduction. Another semi- numerical approach is the one described in . It relies on 
the fact that every one-loop amplitude can be represented in terms of building blocks which 
are one- and two-dimensional parameter integrals in a form which is suitable for numerical 
integration. 

An alternative method is to obtain loop amplitudes by using unitarity to sew together tree 
amplitudes [521 ISHl El] • A difficulty of this approach has been to determine ambiguities of 
rational functions which are present when calculating QCD amplitudes. However, the applica- 
tion of twistor-space inspired methods to one-loop amplitudes led to new insights in this 
context |66j-|70jj, and a rapid development in this direction may be expected in the future. 

Nevertheless, to calculate one-loop amplitudes involving massive particles, as well as to 
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handle infrared divergences due to massless particles, we still have to rely on more conven- 
tional methods. As the size of the expressions for such amplitudes increases factorially with 
the number of legs, efficient methods of tensor reduction become more and more important. 
Although many in principle viable approaches exist, computations which rely on conventional 
reduction methods may get stuck due to the combinatorial growth of intermediate expressions 
for moderate values of N [N ~ 6) already. It is the sheer size of the expressions, together 
with spurious denominators, which in the end hampers a successful, i.e. numerically stable, 
evaluation of the amplitude. Reduction algorithms in momentum space generically lead to 
so-called inverse Gram determinants which vanish if an exceptional kinematical configuration 
is approached. Reduction algorithms in Feynman parameter space in principle overcome this 
problem, but other kinematical determinants are still present in the denominator and one has 
to deal with scalar integrals in higher dimensions |42t 1^ 1321 ISl] • 

In this article, we propose an algorithm which is similar to the one given in [SU] , but improved 
in several respects. First, the new algorithm is designed to restrict and control the occurrence 
of inverse Gram determinants. Second, the formalism is valid for massive as well as massless 
particles, the soft and coUinear divergences being regulated by dimensional regularisation. Our 
method is valid for arbitrary A^, and we give a constructive recipe how to deal with the cases 
where kinematic matrices are not invertible. In addition, we prove explicitly in this formalism 
how A^-point integrals with A^ > 5 in more than n = 4 — 2e dimensions drop out of any physical 
one-loop amplitude. Moreover, we elaborate on the numerical evaluation of the basis integrals. 
Further, the new method is formulated in a manifestly shift invariant way and thus avoids a 
proliferation of terms due to shifts of the loop momentum when the tensor reduction is applied 
iteratively. 

The outline of the paper is as follows. In section |21 we present a non-technical overview 
of our approach, which serves to point out its main features. The following sections contain 
a detailed description of the formalism. In section 01 we define our notation and the general 
setup. The method of tensor reduction by subtraction is described in section |3 In section El we 
elaborate the algebraic evaluation of the building blocks of our reduction. The case A^ = 5 is 
particularly interesting, and in section (HI we give form factors for A^ = 5 which do not contain 
higher dimensional 5-point functions and are free from inverse Gram determinants. The explicit 
proof that these integrals drop out and how the inverse Gram determinants cancel for A^ = 5 
is rather technical and is provided in appendix O In section |71 we deal with the numerical 
evaluation of the basis integrals, by means of multi-dimensional contour deformation, and 
we present explicit checks of the numerical stability near exceptional kinematical situations. 
Section |H1 contains guidelines for the practitioner who is less interested in the mathematical 
details on how to implement the formalism directly into a computer code, before we conclude 
in section ini In the appendices, we provide explicit formulae and useful relations for the direct 
application of our algorithm to multi-leg calculations. 

2 A brief overview of the method 

Before entering into the mathematical details of our formalism we would like to give a short 
overview of the method. 
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We consider one-loop A^-point diagrams with external momenta pi, . . . ,Pn- They are typically 
expressed in terms of integrals in momentum space, with and without loop momenta in the 
numerator. Any algebraic approach to evaluate these diagrams starts by reducing these tensor 
and scalar integrals to simpler objects, with the price to pay that the number of terms increases 
at each reduction step. We will distinguish reduction formulae for tensor and scalar integrals 
in the following. 

Before applying tensor reduction formulae which typically increase the complexity of an expres- 
sion, reducible terms might be cancelled. Numerators of Feynman integrals which contain scalar 
products between loop momenta and external momenta are called reducible if they can be ex- 
pressed by differences of inverse propagators of the given Feynman diagram and by kinematical 
invariants. The momentum representations of the remaining irreducible tensor integrals are 
converted to linear combinations of form factors and Lorentz structures. In our approach this 
is done in a non-standard way, as we express the numerators of the tensor integrals in terms of 
propagator momenta instead of the loop momentum solely. Thus we define a generalised rank 
r tensor integral by 

where Qa = k + Va, and is a combination of external momenta. The method is defined in 
n = 4 — 2e dimensions and thus is applicable to general scattering processes with arbitrary 
propagator masses. Taking integrals of the form as building blocks has two advantages: 1) 
combinations of loop and external momenta appear naturally in Feynman rules, 2) it allows for 
a formulation of the tensor reduction which manifestly maintains the invariance of the integral 
under a shift — > fc + rg in the loop momentum. Such a shift can be absorbed into a redefinition 
of the Tj, Tj — > Tj — tq. The Lorentz structure of the integral (Q) is carried by tensor products 
of the metric g'^'^ and the difference vectors A^- = — r^, which are invariant under such a 
shift. The fact that the sums qa = k + ra, which generically appear in loop diagrams, are not 
split into loop and external momenta, as well as the explicit shift invariance, have the virtue of 
leading to a reduced number of terms in the expressions for the loop graphs. 

A key point of our method is to reduce these tensor integrals by adding and subtracting terms 
such that a reduction into infrared (IR) finite and reduced (or pinched) integrals is achieved. 
By applying the reduction again to the potentially IR divergent reduced terms one generates 
iteratively a separation into IR finite and IR divergent terms. The latter are IR divergent 
3-point functions which can be evaluated analytically in a closed form and separated from 
the finite part of the amplitude. This immediately provides a starting point for a numerical 
evaluation. 

We further show that our reduction of A^-point tensor integrals for A^ > 6 trivially maps to 
5-point tensor integrals. The 5-point case needs special care, as was noted earlier [IHl- We will 
show that all tensor 5-point functions can be reduced to some basis integrals without generating 
higher dimensional 5-point functions nor inverse Gram determinants. The form factors for all 
nontrivial tensor structures for up to rank five 5-point functions will be provided explicitly in 
terms of our basis integrals. 

These basis integrals, i.e. the endpoints of our reduction, are 4-point functions in 6 dimensions 
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/|, which are IR and UV finite, UV divergent 4-point functions in n + 4 dimensions, and various 
3-point functions, some of them with Feynman parameters in the numerator. This provides 
us with a very convenient separation of IR/UV divergences, as the IR poles are exclusively 
contained in the triangle functions. Explicitly, our reduction basis is given by integrals of the 
type 



^jl ■ ■ ■ ^jr 



1 3 3 



iSnM ^ -r (2 - f ) / n - E 

and J^, I^'^'^, ^4^^, 12'^'^ with no Feynman parameters in the numerator. Of course, 2-point 
functions also have to be considered. 

It turns out that for an arbitrary A^-point amplitude, calculated in a gauge^ where the rank 
can only be less or equal to the number of external legs, it is sufficient to consider r < 3 in 
the integrals above. Note that the (n + 2)-dimensional box integrals, being neither IR nor UV 
divergent, can be evaluated using n = 4. The IR divergent three-point integrals are easy to 
handle, and we will give a complete list of all required 3-point integrals for the case of massless 
propagators in appendix iBl 

Our reduction formalism is designed such that any A^-point amplitude can be written as a 
linear combination of the basis integrals without encountering inverse Gram determinants. We 
would like to emphasize that we do not only avoid Gram determinants of rank four matrices 
in the reduction of five-point functions, but obtain form factors which do not have inverse 
Gram determinants from lower rank matrices either. Our form factors are completely free from 
any inverse Gram determinants. In addition, we avoid the proliferation of higher dimensional 
integrals by choosing a convenient set of basic functions, while controlling at the same time 
the occurrence of arbitrary inverse Gram determinants. These two restrictions define our form 
factor representation and lead to the function basis we use. To our best knowledge no explicit 
form factor representations with these special properties have been derived in the literature 
before. To evaluate the basis functions, we propose here two complementary approaches: a 
purely numerical one and an algebraic one. 



^Here we have in mind a gauge fixing for which gauge boson propagators are proportional to g^^. In 
more general covariant gauges where the tensor structure of gauge boson propagators depends on the vector 
boson momentum, the extra gauge dependent terms might lead to higher dimensional integrals with non-trivial 
numerators. However, using the pinch technique |76j . it can be shown that a reorganisation of the integrands 
allows to cancel the gauge dependent contributions at the level of the integrands before the loop momentum 
integration is performed, so that this complication can be avoided in these cases, too. 
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In the algebraic approach our basic buildings blocks, given in eq. 0, are further reduced to 
scalar integrals using recursion formulae. However, this introduces scalar integrals of dimensions 
higher than n+4 Applying scalar reduction formulae, the latter can be remapped to 

(n+2)-dimensional scalar box integrals and n-dimensional scalar two- and three-point functions. 
Doing so, the price to pay is the occurrence of inverse Gram determinants. These dangerous 
denominators are spurious: the expressions can be organised in such a way that inverse Gram 
determinants are multiplied by linear combinations of scalar integrals which also vanish in the 
case of exceptional kinematics, as has been done for example in jT3] . Without such a grouping 
of terms, cancellations between singular pieces in general pose numerical problems, and even 
respecting such a grouping of terms does not guarantee numerical stability. 

If one aims for compact algebraic expressions of loop amplitudes, it is a good strategy to 
express the amplitude in terms of simple scalar integrals and to enforce compensations of Gram 
determinants algebraically. Experience shows that compact expressions can be achieved in this 
way (for examples where algebraic reduction was successfully applied, see jHH IZ21 ESI)- 

In order to simplify the amplitude representation in terms of finite basis functions and to 
enforce explicit cancellations of Gram determinants, it can be useful to exploit relations between 
determinants and sub-determinants (minors) of the kinematic matrix S present in eq. 0. We 
provide a collection of relations useful for this purpose in appendix El Many of these relations 
also served to obtain compact and convenient form factors for the tensor integrals. 

If the amplitude is too complex, the purely algebraic treatment becomes intractable and it is 
advantageous to avoid the introduction of Gram determinants from the start. For this case 
we give a prescription how to evaluate the integrals in eq. (0) without further reduction. The 
IR divergent integrals, being only 3-point integrals in our approach, can easily be handled 
analytically, as they are simple enough to allow for explicit representations of all possible cases. 
For the finite 3-point and the 4-point functions we propose a new numerical method which 
allows for a direct numerical evaluation. By analytic continuation in Feynman parameter space 
and an adequate multidimensional contour deformation we find a numerically stable integral 
representation of the basis functions. The method is described in detail in section 17.11 In 
section 17^ we compare the two approaches of either using eqs. ^ as endpoints of the reduction 
and then proceed numerically, or reducing algebraically until only scalar integrals are reached. 
The comparison shows that near exceptional momentum configurations the purely numerical 
evaluation is stable, whereas the analytic implementation of the basis functions, containing 
inverse Gram determinants, is not. This is true although the terms were grouped such that 
the coefficients of inverse Gram determinants also vanish in the limit of exceptional kinematics. 
We are lead to the conclusion that if the compensation of inverse Gram determinants is not 
possible algebraically, a tensor reduction scheme which avoids them from the start is preferable. 

Thus we propose a method where the form factors are such that their numerical evaluation 
in the kinematically dangerous phase space regions poses no problem. Away from exceptional 
phase space regions, analytical representations can be used safely. These two complementary 
approaches should guarantee a successful evaluation of very complex multi-leg processes. 
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3 Definitions and notation 



3.1 Feynman parameter representations 

We consider general one-loop A^-point graphs as the one shown in Fig. ^ All external momenta 




Figure 1: General A^-point one-loop graph with momentum and propagator labelling. 
Pi are defined as incoming. Momentum conservation implies 

N 

J2p^ = 0. (3) 

i=l 

For future reference we label each propagator qi by the number i as shown in Fig.^ The ordered 
set containing the propagator labels is denoted by S. In Fig. [T]one has 5* = {1, 2, ... , A^}. 

The propagator or internal momenta are labelled accordingly by = k + Vi, where k is the 
momentum running in the loop, and the momenta rj are defined such that Pi = ri ~ rj_i, 
(i = 1, . . . , A^), , To = rjy. Thus one has Qi = Pi + qi-i {qo = Qn)- By momentum conservation, 
one can choose one of the vectors to be zero. Most reduction algorithms specify either r^r or 
Ti to be zero. 

The momentum representation of the scalar A^-point integral in n dimensions is denoted by 

/d"'k 1 

In the following, we will use the shorthand notation dk = d^k/m"^/"^ for the integration measure. 
The ordered set S appearing here as an argument uniquely defines the one-loop integral. We 
will use the set S" as a basic object throughout the paper. 

After having introduced Feynman parameters and performed the momentum integration, I^{S) 
can be written as 



„ N N 

/ I[dz,6{i-J2^i) [RT'"" 

i=l 1=1 



= — - z ■ S ■ z — iS = —- Zi SijZj — iS . (5) 

The kinematic matrix S is defined by 

Sij = {ri -rjf- m\ - m] . (6) 

In general, a one- loop A^-point amplitude will contain A^-point integrals as well as (A^ — 1), (A^ — 
2), . . . , (A^ — M)-point integrals with tree graphs attached to some of the external legs of the 
loop integral. The latter are characterised by the omission of some propagators (say ji, . . . , jm) 
of the "maximal" one loop A^-point graph. They consist of A^ external particles and M < N 
internal lines, where M denotes the number of elements in the set S* \ {ji, . . . , jm}- We give 
some examples in Fig. |21 The corresponding kinematic matrix is denoted by 

= • (7) 

It is obtained from S by replacing the entries of the rows and columns ji, . . . , jm by zero. 




Figure 2: Graphical representation of pinch integrals. Each topology defines an ordered set S. 
The two diagrams correspond to A^ = 8, M = 4, S" = {1,3,6,8} (left), and A^ = 6, M = 4, 
^ = {1,3,5,6} (right). 



In this way one can keep track of the pinching of propagators in the iterative application of 
reduction formulae without changing the labels of the rows and columns of reduced matrices 
g{h --im} -^yith respect to the maximal set S. 



For example, for 5* 
/ 

S = 

V 



{1, 2, 3, 4} one has, with Sj = p| and Sij = {pi + pj] 



-2 ml 



S2 



S2 ■ 

si ■ 



2 



rrir 



2 2 



2 



nil " ^2 
-2ml 

S3 — 7712 ~ ""^3 



m^ 



m. 



S23 

53 ■ 

54 - 



- mi 

2 

-2ml 

2 

mi 



-ml 
ml 



m 



S12 - m^ 

Si 



ml 
- ml 



mi 



2 2 
- ?Ti3 — m| 

-2 ml 



J 



S12 - "^2 

The symmetric (4 x 4) matrix S^'^'^\ which corresponds to the pinching of propagators 
4, is now defined by 

/ -2 ml 




^{2,4} 



S23 - 
\ 



mi 










■523 



2 2 

m^ — mg 





-2 ml 









2 and 



(9) 
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Inverse matrices are labelled analogously. Although S^^^'"^"^^ in our definition is not a regular 
matrix it has a uniquely defined pseudo-inverse. We recall that the so-called Moore-Penrose 
generalised inverse P to a symmetric matrix S is uniquely defined by the properties jTTl [TH] 



vsv = V , svs = s , vs = sr. 



(10) 



This concept will also be used below. To construct the pseudo inverse here one simply has to 
invert the sub-matrix of S^^^'"^"'^ with the zero rows and columns omitted and promote the 
result back to an x iV matrix by inserting zeros for the rows and columns {ji, . . . , j'^}. In 
our example one finds, with A(x, y, z) = + y"^ + ~ 2xy — 2yz — 2xz: 



C{2,4}-1 _ ^ U U U U 





/ 


2 ml 





S23 - ml - ml 


\ 


1 
















\{s23,ml,ml) 




2 2 
— mj — mg 





2 ml 







\ 











0/ 



In the following, S ^{S \ {ji, . . . ,jm}) = S^^'^'''''^"'^ ^ has to be understood in this sense. 

Using these conventions, Feynman parameter integrals with propagator pinches and Feynman 
parameters zi-^ . . . zi^ in the numerators can be defined as 



1) 



Jr] S \ {ji, . . .,jm}) 



^ N N 
i=l 1=1 



2\n/2-N 



(12) 



Whenever one index of the pinch set Q = {ji, . . . , jVrt} coincides with one numerator index in 
the set JV = {li, ... Jr} this integral is trivially zero: 



If MnQ^U 



I^i^^;S\Q) = o . 



(13) 



The above conventions lead to simple expressions in our formalism and are designed for iteration 
purposes, as will become clear below. 



3.2 Definition of form factors 

Algebraic expressions of an amplitude typically consist of spinors and Dirac chains of momenta 
and polarisation vectors, depending on the process under consideration. The loop momenta 
which appear in the numerator of a Feynman diagram often come in the combination = k+Vi, 
like for instance in fermion propagators. We keep this natural combination in our tensor 
reduction formalism. Therefore we define tensor integrals as 

By setting ai, . . . , = A^, and using momentum conservation to set rjv = 0, we can always 
retrieve the commonly used form 

rr f 2"' 2^;, ■ (15) 
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In this more "conventional" approach, one of the gf is specified to be k^, which singles out one 
propagator and defines a standard form. After a reduction step one obtains integrals which 
are not of standard type, such that a shift operation A; ^ A; + is necessary to remap to 
the standard form, giving rise to 2*" terms for a rank r tensor integral. Our formulation, being 
manifestly translation invariant and thus making such shifts obsolete, avoids here a proliferation 
of terms. 

As pointed out above, shifts of the loop momentum can be absorbed into a redefinition of the 
Tj vectors. To achieve a manifestly translation invariant formulation, we need vectors which 
are invariant under shifts r|" + r^^. This motivates the definition of the shift-invariant 

vector A^-: 

Af, = rf - r>; = gf - (16) 

Apart from metric tensors g'^'^, the Lorentz structure of the integrals will be carried by these 
vectors. 

As will become clear below, we have to distinguish the cases < 5 and A^ > 6. In the case 
A^ < 5 we will express the different tensor structures in terms of metric tensors and difference 
vectors, A^. Tensor integrals are expressible by linear combinations of such Lorentz tensors 
and form factors denoted by A^^l^^^ (S), bI^'^{S), C^,''^{S). A^'^ is the coefficient of the Lorentz 
structure containing only difference vectors. B^'"^ belongs to exactly one metric tensor and 
(r — 2) A^- vectors, and C^''^ is the coefficient of the Lorentz structure containing products of 
two metric tensors. Thus our form factors for A^ < 5 are defined by the formula 

jl---jr-2eS 
jl--jr-4eS 

where the distribution of the r Lorentz indices /ij, and momentum labels to the vectors A^'^., 
denoted by [■ • in eq. (|T7jl . is illustrated in the following equations (|T!?jl to ((221) ■ 

lUS) = (18) 



ir'{a,;S) = J2^f'a.A^'\S) (19) 
/™(ai,a2;5)= Ka. K.. KH^) + 9'' B'^'^S) (20) 



ir-'^'^iau a2, a,; S) = K\ Ka, K'as KLi^) 

+ J2 ir'^ ^fL + 9'"-' Afa, + 9''"' A^) B^'\S) (21) 
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I^^^^'^-^^^ia,, a„ as, a,; S) = Ka. Ka. Ka, KLiS^) 



_|_^M2M4 A/'l A/^3 _|_ „M3M4 A A*! AM2 \ D^'^/rtN 



fl4 



_|_ (^glJ.llJ.2 ^M3M4 _j_ ^/ilM3 ^^2/^4 _j_ ^/i2/i3 gfJ-lfJ-i^ C'^''^{S) (22) 

We recall that standard form factor representations can be simply obtained by replacing aj = N 
for all j, together with r^r = 0. This also shows that the form factors do not depend on the 
introduction of the difference vector. The form factors are shift invariant by themselves. One 
main result of the paper will be the explicit representation of all these form factors in terms of 
higher dimensional 4-point parameter integrals and n = 4 — 2e dimensional 3-point parameter 
integrals with nontrivial numerators, see eq. (P)). 

For > 6, the tensor reduction will be done in such a way that only the form factors for 
= 5 appear. The Lorentz structure of A^-point rank r tensor integrals does not require 
the introduction of additional factors of g'^'^ as compared to the N = 5 case, only additional 
difference vectors A^- appear. This is due to the fact that for > 5, four linearly independent 
external vectors form a basis of Minkowski space. We note that for A^ = 5, one could already 
express the metric by external momenta, but this would introduce inverse Gram determinants. 

Before closing this section we would like to note that momentum integrals, Feynman parameter 
integrals and form factors are naturally related by | ^ I44 | IHH] : 



/---(a,,..., a.; 5) = {-lY ^2 E liar- ^n- ' ' ' ^^r]^'^' 

m=0 ^ ^ jl--jr-2m = l 

Xr/'"^{j,...,Jr.2m;S) (23) 



where ■ ■ ■ ,jr-2m ; S) is defined in eq. (fT^ . In eq. (j^ . [r/2] stands for the nearest 

integer less or equal to r/2 and the symbol ®m indicates that m powers of the metric tensor 
are present. It is obvious from this formula that the ^f'^^'-terms are always associated to integrals 
in more than n dimensions. 

In the following sections 0] and we will formulate a reduction formalism for tensor integrals 
of the type (fT^ and give convenient representations for the form factors defined in eq. ()17|1 . 
In section O we will show that higher than n = 4 — 2e dimensional integrals for A^ > 5 can 
be avoided completely if the external kinematics is defined in 4-dimensions, a fact which was 
already carefully investigated elsewhere |121 HEl HO] • 



4 Tensor reduction by subtraction 

The formalism described in this section naturally leads to a separation of IR divergent and 
finite expressions and does not produce spurious Gram determinants. The reduction is based 
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on a subtraction technique which is analogous to the one used in ^Dj for the scalar case. Before 
we come to the tensorial case, let us recall the procedure of jSH] for the scalar case, recast into 
the notation of this article. 



4.1 Subtraction for scalar integrals 

Our aim is to split a scalar A^-point integral as defined in eq. into an IR finite part and a 
possibly IR divergent, but simpler part. Therefore we make the ansatz 



(24) 



We can see that Idiv is a sum of reduced integrals where one propagator has been pinched. 
Now let us consider I fin after having introduced Feynman parameters. In order to arrive at a 
quadratic form in the loop momentum, we perform the shift 

k = l-^ZiTi. (25) 

i&S 

This shift transforms the denominator to the form — R^, where is defined in eq. (0). For 
the numerator, we use 

■ Aki = 2^Sii + Sjk — Sik — Sji) , 



to obtain 



l-Y^.{S){qj-mi) = -(/2 + /?2) 5^6,(5) + 5^2, 1-5^6,(5) {5,, +2/ -A,,} 
ies i£S jes L 



(27) 



The term linear in the loop momentum / vanishes due to symmetric integration and we conclude 
that if the equation 

YHS)S^, = 1 , j = l,...,iV (28) 

ieS 

is fulfilled, the term in square brackets in eq. (|77|) vanishes. Then I fin is given by 



/•I r d"^ 

IMS) = -B{s)r{N) / l[dzj{i-J2'^) 



d^^l P + R^ 



V2 {P-R^)N 

-BiS)iN-n-l)I^-^'\S) (29) 



es les 

n+2/ 



B{S) = 5^6,(5). (30) 



iG5 
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If it is clear from the context which set 5* we refer to, the argument (5*) is omitted in B{S) 
and bi{S). If the bi belong to a reduced kinematic matrix S^^^ where the j^^ row and column is 
zero, associated with the set \ {j}, one has bi(S\ {j}) = J2kes\{j}('^^''^)ki ■ -^^^ simphcity of 
notation, we introduce the shorthand bi{S\ {j}) = and correspondingly B^^^ is defined as 

B^^^= J2 b?^- (31) 

i&S\{j} 

All that remains to be shown now is that eq. ()28|1 indeed has a solution for the reduction 
coefficients bi for arbitrary A^. 

In the case of non-exceptional 4-dimensional kinematics^, rank{S) = min{N, 6) holds. Thus S 
is invertible for < 6, and eq. (j28|) has the unique solution 

h = Y.Sk'- (32) 

fce5 

If S is not invertible, we proceed as follows. First we single out the a*^ row and column in S, 
where a is an arbitrary element of the set 5, to write Sij as: 

5,, = -Gf^ + Vt'^ + y/") (33) 

with 

g!;) = . A,., V^('^) = Ai-m^ (34) 

The Gram matrix G^-"'' is understood as an A^ x A^ matrix. By definition its entries are zero 
whenever i = a ot j = a. For ease of notation, we will omit the superscript (a) in the following 
if it is clear from the context. Using the above relations and distinguishing the cases j = a and 
j 7^ a, eq. (PHj) is equivalent to the two equations 

b,Gi, = B{V,-Va) (35) 

i€S\{a} 

J2 b^ = l-2BVa. (36) 

i£S\{a} 

Eqs. fj35p .()36 p may be solved in the following way. As G is not a regular matrix we first construct 
the generalised inverse H of G, defined as in eq. fjlOj) . To this end we introduce four linearly 
independent 4- vectors Ej^^^ ... ^ forming a basis of the physical Minkowski space, and a 4 x A^ 
coefficient matrix TZ, such that 

4 

Afa = E ^™ = ■ • (37) 

m=l 

^We call a kinematic configuration defined by N external momentum vectors exceptional if the Gram deter- 
minant built from a subset of four of these vectors vanishes. 
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Here T^ma = for all m. The Gram matrix G is expressible as^ G = TiT^GlZ. The matrix IZ is 
of rank 4, and thus the 4x4 matrix TZIZ^ is invertible. The sought matrix H is then uniquely 
defined by 

H = n^{nn^)-^G-\nn^)-^n . (38) 

Before going on we note that the case of exceptional kinematics is now easily dealt with. If 
the external vectors only form a space of dimension k less then four, one may choose k basis 
vectors i?/^^ ... defining 7^ to be a x matrix. This is the only change to be made. 

The system of equations admits solutions if and only if the consistency condition 

B {In-i~GH),,6vI^^ = (39) 
jes\{a} 

is fulfilled, where we have defined 

^ jeS\{a}. (40) 
We will not always write the sum J2j(^s\{a} explicitly in the following, but denote products like 

^,^S\{a}K^,V,hJ{K.V).. 

Since, for > 5, in general GH ■ 6V ^ bV ^ a solution of ()39|) for > 5 exists if and only if 
5 = 0. The solution spans an (iV — 5)-dimensional space which is just the kernel of the Gram 
matrix. It can be parametrized by (A^ — 5) vectors [/(^'---.a^-s)^ ^g^- ^ote that 

GH = HG = n^inn^Y^TZ (41) 

and that the projector onto Ker(G') is given by 

K = -7^^(7^7^^)-^7^ . (42) 

It follows from the definition of H that K ■ SV G Ker(G'). Now one can choose [/(^~^) = 
K-5V/{6V-K-5V) parallel to 5V and the other [/-vectors orthogonal, 5^^. f/(fc=i.-.^-6) = o. A 
general vector in Ker(G) is then parametrised hj U = Ylk=i f^kU^^'' + af/^^"^-' and a solution 
to eqs. is given by 

6a = - 5^ &A:. (44) 

keS\{a} 

We recall that this solution is valid for all A^ > 6. For exceptional kinematics and A^ < 6 an 
analogous solution can be derived. In the other cases one can simply use eq. ()32|1 . 

•^Note that the pseudo-inverse of G and the singular value decomposition of G given in ref. [23 are related 
in a simple wav |78|: Defining an orthogonal 4x4 matrix O which diagonalises Gmn (such a matrix O al- 
ways exists for a symmetric matrix G), i.e. {O^GO)mn = i^mSmn, one has Gij = {TZ^ O G O TZ)ij = 
(^"^ 0)ik ijJk [O^ T^)kj = Sfc=i '^ik vj^j, which means that the Uik and v^^ used in ref. are related 
to the matrix TZ given here by Uik — {TZ^O)ik and vj^- = {O^TZ)kj- 
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Near some special momentum configurations the denominators of the reduction coefficients, 
6V ■ K ■ 5V or det(5), may cause numerical instabilities in realistic applications. In contrast to 
inverse Gram determinants which can be viewed as technical artefacts of the tensor reduction 
down to scalar integrals in momentum space, the singular behaviour of the latter is due to 
physical singularities like soft /coUinear configurations or thresholds and therefore they are 
unavoidable. We would like to point out that this is not a problem related to the use of the 
pseudo-inverse. If one would use the singular value decomposition, this problem would occur 
as well. Experimental cuts very often exclude these numerically dangerous phase space regions. 

Finally, we quote the interesting relation between the kinematic determinants 

5 det 5 = (-1)^+^ det G^'^) , (45) 

N. 



where G^^ is the Gram matrix defined above. The relation is valid for all a G and arbitrary 



We conclude that in the case > 6 one always has -B = 0, if the external kinematics is defined 
in four dimensions. This means that Ifin{S) = in eq. ()29|) . In the case N = 5 one finds 
IfiniS) = B{n-A) I^^\S). As I^^^S) is IR and UV finite, the whole term is of 0{e) and 
can be dropped in phenomenological applications. For = 4 we see that 4-point functions can 
be represented in terms of 3-point functions and (n + 2)-dimensional 4-point functions. As the 
latter are IR finite, a splitting into IR divergent and IR finite integrals is achieved. We would 
like to emphasise that by iteration of eqs. and (j^^ an arbitrary scalar A^-point function 
can be algebraically reduced to 3-point functions and (n + 2)-dimensional 4-point functions. In 
this representation IR divergent integrals are naturally separated from finite contributions. 

As an example we give here the 6-point function in terms of (n+2)-dimensional 4-point functions 
and 3-point functions. In our notation, up to 0(e) terms: 



jes k€S\{j} 



{j} 

k 



B^^'>'Ur'{S\{j,k})+ E bj'''hUS\{j,k,l}) 

ies\{j,k} 



From this representation, IR divergences can be trivially isolated. 



4.2 Subtraction for tensor integrals 

We now extend the above reasoning to the tensorial case, i.e. we will split a general tensor 
integral into an infrared finite part and a part which contains possible IR divergences, the latter 
having one rank less and one propagator less. We will give a constructive algorithm how to 
proceed for any number A^ of external legs. 

Let us write eq. (HH) as follows: 



jn, 111... fir 



(ai, ...,ar;S) 



dk 



E 

i65 



dk 



Mar 



(46) 
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The last line corresponds to (A^ — l)-point tensor integrals of rank (r — 1). The coefficients C^^^ 
will be determined such that the ffist term in eq. ()46|) is an IR finite expression. 



As in the scalar case, one introduces N Feynman parameters Zi and makes the substitution 
()25j) to obtain the form — for the denominator. Under the shift (j2SI), the momenta qa 
become 

ies 

Now let us consider the vector A'^ , given by the square bracket in the first line of eq. ()46|) 



ie5 



2 2\ 



(48) 



where Cj^ is defined by the following equation, which will be the cornerstone of our derivation 
for general N : 

5^5,,C;„ = Af„, aeS. (49) 

Note that this equation is analogous to eq. (PH|) for the scalar case. We will solve equation 
in the next subsection for arbitrary kinematics and arbitrary A^. 



The vector A'^ transforms under the shift as: 



i2 2 



.J6S 



(50) 



where eq. (pUj) and the definition 



iG5 



k a 



(51) 



k€S 



have been used. Thus we see that if eq. fj49p is fulfilled, all the terms contained in A'^ are either 
proportional to the loop momentum / or to R^. 

Provided that eq. (jiUI) is fulfilled, and using Aij = Aia^ + Aa^j, where 02 G S" is arbitrary, 
eq. ()5Up can also be written in the following form : 



V ies J 



where 



ai 02 



J a.2 



(52) 



(53) 
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From the expression (j52p one can see immediately that the first term on the right-hand side 
in eq. ()46|) has no infrared divergences: After the shift ()25|) . the integral in the first line is 
proportional to / or + and, after integration over the loop momentum, will give some 
higher-dimensional integrals (with Feynman parameters in the numerator for r > 1). The term 
in the second line of eq. ()4(jj) is divergent but has one tensor rank less and one propagator less. 
Thus, if eq. fHnjl is fulfilled, the ansatz leads to the desired splitting of a rank r A^-point 
integral into IR finite terms and (iV — l)-point integrals of rank r — 1. 

4.3 Solving the defining equation for arbitrary 

It remains to be shown that eq. (j49p has a solution for general A^, and to construct such a 
solution explicitly. 

As in the scalar case, the solution is simple if S is invertible, i.e. in the case of non-exceptional 
kinematics for A^ < 6. In this case, eq. has the unique solution 

kes 

On the other hand, if A^ > 7 or in the case of exceptional kinematics, S is not invertible, so 
eq. (jinjl does not have a unique solution. However, an explicit solution can be constructed in 
the same way as has been done in section l^?T] for the scalar case. To this end, we again write 
the matrix S as in eqs. Using the definition (|in|) for 6Vj°'\ eq. may be rewritten 

as 

ies 

G^Ct, = -A';^ + 5V}^^VI:, (55) 

ieS\{a} 

J2 = A;:,-21^Wv,^ (56) 

i€S\{a} 

where also '^i^s^ib ~ "^b been used. Eqs. (jHSj) , (jSSl) can be solved using the same pseudo- 
inverse as already constructed in section I01 The system of equations admits solutions 
if and only if the consistency condition 

J2 i^N-i - GH),, (a;, - SVj''^ V^) = (57) 
jes\{a} 

is fulfilled. Since J2jes\{a}i^N-i ~GH)ijAj^ = whereas, for A^ > 5, in general GH -dV ^ 5V , 
a solution of for A^ > 5 exists if and only if = 0. The general solution of for A^ > 5 
is thus given by 

Ctb = - Yl H.,A';^ + Wt , zeS\{a} (58) 
jes\{a} 

'■^ab — 2^ ^jb ) ^b ■ 

iGS\{a} 
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The vectors Wj^ span the (A^ — 5)-dimensional vector space Ker(G), the kernel of G. To 
construct a basis of Ker(G) we again use the vectors introduced in section I^^Tl The 

vectors Wj^ are then parametrised as 

Wt = f3'M-,{K-6V), + J^(3i'uP . (59) 

1=1 

Substituting the parametrisation into eq. (j^Ujl yields: 

^N-. = 6V.K.6V (^"^+ ^ SV.H^.aA (60) 

\ ^d&S\{a} ) 

... Ar-6 — arbitrary 4-vectors . 

If = 6, G is not invertible, yet the solution is still unique, as it has to be since S 

is invertible. If > 7, neither G nor S are invertible. In this case ()49|) still admits solutions, 
which are however no more unique, but span the [N — 6)-dimensional affine space defined by 
(jKH|) - (jHIH) . It also has to be emphasised that the above construction is equally valid if the 
external momenta become linearly dependent not due to the fact that > 6, but due to an 
exceptional kinematic configuration. 

With relations (|37|) and (|58|) at hand, it is now easy to see that for A^ > 6 one has 

^C;,A^^, = -^gi-^ and thus T^^C = . (61) 

The subscript "[4]" in indicates that T^^ is 4-dimensional and not n-dimensional here. 

Further, we already saw that 

V^^ = for A^ > 6 . (62) 

Relations ()6H) and ()62|) have an important consequence. They mean that A'^ in eq. ()52j) is zero 
and thus they imply that no higher than ri-dimensional A^-point functions with A^ > 6 can be 
generated by the reduction: In eq. ()46p . only the pinched terms survive, leading to 

Jr^-^^(ai,...,a.;^) = - C^^^ j;'^r"^(a2, . . . , a.; ^ \ {j}) (AT > 6) . (63) 

In this sense the tensor reduction of A^-point integrals with A^ > 6 is trivial: Integrals with 
A^ > 6 can be reduced iteratively to 5-point integrals, without generating higher dimensional 
remainders. Therefore form factors for A^ > 5 are not needed. 

The case N = 5 needs special attention and was already carefully analysed in the Feynman 
parameter approach ^H]- We will give a full analysis of the 5-point case within our formalism 
in section lO 

In order to make contact to ref. |50j, we note that the reduction is equivalent to applying first 
eq. (32) and then eq. (B.7) of [20], but the method advocated here gives very naturally certain 
algebraic relations between reduction coefficients, which have been exploited in order to avoid 
Gram determinants (i.e. factors of 1/i?) as far as possible. 
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5 Form factors for = 3, 4 and algebraic representation 
of basis integrals 

In the following we will apply the formalism explained above to provide explicit formulae for all 
form factors present in the reduction of 3- and 4-point integrals. The form factors for N = 5 will 
be given in section El and the 2-point functions will be listed in appendix ^ The form factors 
will contain scalar integrals with non-trivial numerators. In this representation no inverse Gram 
determinants are present, which serves as an ideal starting point for a numerical approach^. 

Nonetheless we will also provide an analytical representation for the basis integrals in terms of 
scalar integrals with trivial numerators only. The price to pay in this case are inverse Gram 
determinants. For the calculation of a cross section this representation is numerically stable 
in most parts of the phase space. However, numerical problems can arise at the kinematical 
boundaries. In section [7| we will propose a numerical evaluation of the basis integrals which 
also works for exceptional kinematical configurations. 

5.1 Three-point integrals 

The form factors for the 3-point tensor integrals, following directly from eqs. (fT7j) and ()23|). are 
given by 

Af\S) = -m;S) 

KLiS) = -ink,k,h;S) (64) 

If one would like to express the integrals with nontrivial numerators above in terms of scalar 
integrals only, the formulae given below in eqs. to (fTO)) can be used if S is regular. The 
singular cases occur if either the 3-point function is IR divergent or one hits an anomalous 
threshold. Anomalous thresholds appear in scattering processes with unstable external par- 
ticles. In very special kinematic situations all internal particles in a given diagram can go 
on-shell, which can cause integrable singularities inside the phase space [HDl- One can show 
that in the 3-point case det(iS) = is a necessary condition for an anomalous threshold. For 
the case where IR divergences are present, we give a complete list of 3-point integrals with 
massless propagators in appendix iBl 



^3(^1;^) = I 



l3iS)-J2b,I^{S\{j}) 

jes 



+ E'^^u^2(^\{j}) (65) 

jes 

^Parameter integrals with non-trivial numerators, corresponding to tensor integrals, are also dealt with 
in j7H| . where n-dimensional tensor 3-point and 4-point functions with massless propagators are expressed by 
hypergeometric series expansions. 
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Isih, h, h] S) 
ir\h;S) 



-sfjjr'is) + hM - s) +J2^rj m^; s \ {j}) m 

jes 



jes 
1 1 
B{n-2 

1 



'ms)-J2bims\{i}) 



(67) 
(68) 



les 



I2; S) 



jes 
jes 

i- I^+'iW, S) + I^+'{h; S) + I^h, W, S) 



-J2bjI-{hj2;S\{j}) 

jes 



(69) 



(70) 



The two-point functions are given in appendix 1X1 The scalar three-point function I^{S) is 
well known, see for example |^ 1741 175] . 

By iterating the above formulae one arrives at a representation in terms of scalar integrals with 
trivial numerator. As a result of the iteration, inverse Gram determinants up to the third power 
occur. To improve the numerical stability, the bracketing of the terms as given by eqs. 
and (jUHj) to (f7n|) should be respected. 

5.2 Four-point integrals 

All form factors of the 4-point tensor integrals are expressed in terms of {n + 2)- and {n + 4)- 
dimensional scalar box integrals and n- and (n -|- 2)-dimensional triangle integrals, with up to 
three Feynman parameters in the numerator. 



At'\S) 



Bir\s) + j2bjms\{j}) 

jes 

-biir'is)-j2siims\{j}) 

jes 

+ [^iil ^3 (^1; s \ {j}) + s.,i mw, s \ {j})] 

jes 



(71) 
(72) 

(73) 
(74) 



Bf'\S) 



(75) 



21 



- [bi, /r'(^2, h; S) + bi, I2+\h, h; S) + bi, I2^\k, h; S)] 
-\ E i^ih h; S \ {j}) + Srl^ I^ih, h; S \ {j}) 



3 

jes 

+ Sr^ll-{h,k;S\{j})] (76) 

C^-\S) = ^I2^\S) (77) 

<t(^) = -ll2^\hj2;S) (78) 
^hiisiS^) = /'''(^i' ^2; ^3, ^4) + f''\h, h; k, k) + f^\k. k; k, k) 

+f^''^ik, I3', k, k) + f^'^ik, k', k, k) + f^'^^ik, k', k, k) 

+g'^\k;k,k,k) + g''\k;k,k,k) 
+9''\k;k,k,k) + 9''\k;kj3Ji) 

f^''^{k-,k'-,k-,k) = —-;^Si^]^I2'^'^{k-,k'-,S) 

9''\k;k,k,k) = 6;,/r'(/2,/3,/4;^) + ^ J]5-,;/3"(/2,/3,/4;^\{j}) (79) 

There are six different combinations of terms /^'^(. . .) and four different combinations of terms 
5f^'^(. . .). The combinations are imposed by the symmetry of these objects and represent all 
different distinguishable index distributions. The object /^''^ is symmetric in the first and last 
two indices. The symmetry in the first two indices is manifest, as iS;"}^ is symmetric. The 
symmetry in the last two indices is only induced when the form factors are combined with 
formulae (jl7|) . because the summations symmetrize these indices. The object (7^'^ is symmetric 
in the last three indices in the same sense. 

The form factors given above are free from any inverse Gram determinant and contain only 
n-dimensional 3-point integrals with up to three Feynman parameters in the numerator, {n + 2)- 
dimensional 3-point integrals with maximally one Feynman parameter in the numerator, and 
(n + 2)- and (n + 4)-dimensional 4-point functions with up to three respectively one Feynman 
parameter in the numerator. These integrals form our basis set, i.e. the endpoints of the 
reduction. We will see that for > 4, no additional basis integrals are required. This basis 
set, being free from inverse Gram determinants, serves as an ideal starting point for a numerical 
evaluation. 

Of course, further algebraic reduction down to integrals with no Feynman parameters in the 
numerator is also possible. The following formulae and the results for the 3-point integrals 
given above allow to achieve an algebraic representation of the form factors in terms of scalar 
integrals, i.e. integrals with no Feynman parameters in the numerator, only. Again, such 
a representation produces 1/B terms, which is equivalent to the presence of inverse Gram 
determinants. Exploiting the relations given in appendix \D.1\ the following results can be 
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derived: 



\ J2 ^iims \ {j}) -Ij^^j mi; s \ {j}) I (80) 

jes jes ^ 



I2+\hMS) = ^{hiJl^\W,S) + hiJl^\h-S) (81) 



35 



+ 1 E ^3 (^2; s \ {j}) mh, h; s \ {j}) I 

^ jes ^ jes J 

-I {^iX ^r'(^3; s) + ir\h, s) + srx s)) 

+ \ (E ^^il ^3 (^1' S \ {i}) + E ISih, h; S \ {^}) 

ies ies 

+ E ^3 (^2, h; S \ {^})) - ^ E ^3 (^1, h, h, S \ {^}) I 

ieS ies ^ 



/rV) = _\) ^ { I2^\S) - gfe, ir\S \ {j})| (83) 

/r^(/;^) = -^\kir\s)+ir\i;s)-j2b,ir'ihs\{j})\ (84) 

Applying these formulae and eq. (jHH), an algebraic representation of the form factors for 4-point 
tensor integrals in terms of 3-point functions and six- dimensional integrals without Feynman 
parameters in the numerator is achieved. The six- dimensional scalar box functions for massless 
internal lines are well known and can be found for example in jlUl [^. In section [7| the 
behaviour of the form factors for 5 — in this algebraic representation will be compared to 
the one given by eqs. (f?^ to (f?^ . where inverse Gram determinants have been avoided. 

6 Form factors for N = 5 

In this section we will give numerically stable form factor representation for five point functions. 
We have already seen that the reduction of tensor integrals for > 6 does not produce higher 
dimensional remainders, see eq. ()63|) . The case = 5 is, in a sense, the most complicated one. 
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The derivation is rather technical and can be found in appendixO For five-point integrals up to 
rank three, the absence of higher dimensional five-point integrals already has been demonstrated 
by explicit calculation in {351 • appendix O we show that these integrals drop out for 5- 
point integrals of arbitrary rank and we derive representations which are free from both higher 
dimensional five-point integrals and inverse Gram determinants at the same time. To the best 
of our knowledge, such a representation has not been given in the literature before. 

By application of eq. ()(I103j) and reduction of ?7,-dimensional box integrals 12 to 3-point 
functions and higher dimensional box integrals, we obtain, dropping 0{e) terms: 

j2b,B^^ur'{s\{j}) + j2 E b,bi^h2{s\{j,k}) 

je5 jes kes\{j} 

jes jes kes\{j} 



(85) 

(86) 

(87) 



jes 



( ^ii] K + s-l bi, - 2 SfX h + h I.) i2^\s \ {j}) 



jes 



jes k&s\{j} 



jes 



(88) 
(89) 



4^3 (^) = ttE [^^'(^3;5\{j}) 



jes 



X ( 2 SfX b, - Sr,] - Sjl^ bi, - b, S^^y-\, + l, ^ l, + l, ^ l,] 



jes 

jes kes\{j} 
-f Zi ^ /2 + ^1 ^ h] 



(90) 



C^'^(^) 



n - A 



4 V 3 



jes 

rn+A 



1 + 



E^:'-^^'(^\{^'}) 

jes 



X I2+\S\{j})-b,I2+\h,h;S\{j}) 
i 5r^; S \ {j}) - I Sjl I2^\h- S \ {j}) 



kes\{j} 



(91) 



(92) 
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f^'^^ih, h] hi h) + f^'^ih, h] h, h) + /^'^(^i, h] h, h 



+/^''^(^2, ^3; ^1, ^4) + f^'^^ih, h', h, h) + f^'^{h^ h', h, h) 



(93) 



f^'^ih, h] h, h) 



-2 /r'(/i, k; S \ {j}) ( 2 h, - S^l bu - Sjl k, - h, Si^y-\ u) 
+ 1 H mh,k;S\{j,k}) {S7l^Si^}-\u+S-ilSi^^-\i,) 



g ' (^1; k, k) 



k&S\{j} 



(94) 



The combinations of /^'^(. . .) and g^''^{. . .) appearing in A^^^^i^i^{S) are imposed by the symmetry 
of these objects and represent all different distinguishable index distributions. The object /^'^ 
is symmetric in the first and last two indices and the object g^'"^ is symmetric in the last three 
indices when combined with formulae (fTTjl . 



n- A 



b,I2^\l;S\{j})--Srll2^\S\{j}) 



n 



ir'ih;s\m 



"its 
eE 

X ( 57^; bi, + 57^; bi, - 2 5,;J3 b, + b, S^^^-\, I,) +1,^1^ + 1,^1,^ 

- ^/r'(5\ {j}) {sji^si^^-\,,,+ST^s^^^-\,,,+sj^s^^^-\,,:) 

+ i2^\kMM\S\{j})b, 

+ |i /r'(/i, k; S \ {3}) Sjl + 1,^1^ + 1,^1,^ 

+ \ H [ir\h:s\{j,k]) 



kes\{j} 



X 



(95) 



iE 



{f^'^{k, k, k', k, k) + 9 combinations) 



+ (^'^'^(^i; k] k, k, k) + 9 combinations) 



(96) 
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f'\k, k, h; h, h) = -\ ^3 (^1, h; S \ {j, k}) (57/^ Si^^-\i, + Sj^^ S^^^-\i, ) 

kes\{3} 

-2 /2, /s; S \ {j}) ( Sjl^ + Sr^ b,, - 2 S^X + '^^^'^^'u 

r75'5(/i,/2;/3,^4,/5) = ^/r'(/i,/2;5\{j}) 

X (57^; 5i^i-^3U + Sjl S^^}-\i, + Srl^ S^^}~\i, ) (97) 

The ten different combinations of the functions /^'^(. . .) and g^'^{. . .) are all different distin- 
guishable index distributions. The function /^'^ is symmetric in the first three and last two 
arguments and the function g^'^ is symmetric in the first two and last three arguments when 
combined with formulae (fTTj) . 



7 Numerical evaluation of the basis integrals 

In this section we present a method to evaluate the basis integrals of our reduction formalism 
numerically. The method is an alternative to an approach proposed earlier [HI] , which would 
also be viable, but which needs more analytical input. The method presented here allows to 
deal easily with the case of complex masses, which is necessary if unstable particles are present 
in the loop. 

First we will explain the mathematical details of the contour deformation. Then we will show 
a comparison between the numerical and the algebraic implementation of some basis integrals. 



7.1 Contour deformation of parameter integrals 

The method which will be presented here is based on contour deformation in Feynman pa- 
rameter space. As explained above, in our formalism it is sufficient to evaluate the following 
functions, which are the endpoints of our reduction 

"^0 i=l 1=1 



= -r(e) / \{dzj{i-Y. 

-^0 1=1 1=1 



'~z-S- z/2 - id) 

1 3 3 



Zj^ 



z-S-z/2-i6Y 

1 4 4 



X 4 4 

ir\ji,---,Jr) = r(e) / lldz,6{l-J2^ih — ^ 

-^0 i=i 1=1 [-^-^ 



z/2-i5y 



(98) 



together with integrals of the same type, but with no Feynman parameters in the numerator, 
and two-point functions. We will find numerically stable representations of these integrals as 
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special cases of a completely general derivation which is valid for scalar integrals of the form 

^ „1 iV N 

J A 1 J 1 



(99) 



in the case when no IR/UV divergences are present. The IR/UV singular cases are discussed 
below. 

For loop calculations with unstable particles, it is necessary that internal propagators can be 
defined with complex masses, TVI^ = M| — iMjVj. The denominators of the integrands in 
eq. (PU]) are defined accordingly by changing the quadratic form to 



1 <c 

z ■ o ■ z — I 

2 



^^-\z.S-z-i(^z^ (^z,Mf^ (100) 



The finite width improves the stability of a numerical integration but is not sufficient to guaran- 
tee a stable evaluation. We construct a contour deformation in those parameter integrals which 
belong to propagators with zero width only. This has technical reasons which will become clear 
below. 

In a first step we eliminate the delta function in eq. (jU^ by decomposing the integration region 
into iV sectors where in each sector, one of the Feynman parameters is larger than all the others. 

1 = ^ Q{zi > ^1, . . . , Zi-x, Zi+i, ...,zn) (101) 

1=1 

Note that the splitting into sectors can be avoided by a clever choice of Feynman parametrisa- 
tion adapted to the special case at hand. However, as we are interested in a general, automated 
approach, we do not pursue this option. The basic integral decays into sector integrals J; 

N 

InUu •••,>) = (-l)^r(iV - D/2) MN, D,j,, (102) 

/=i 

The latter can be written as integrals over the (A^ — l)-dimensional unit hypercube. Focusing 
on one sector, say sector I, and introducing the vector T = (ti, . . . , l,ti, . . . , tN-i) which 
defines the — 1 coordinates ti, . . . ,t^_i of the unit hypercube, one gets: 

'^-1 ^ N-D-r T T 

J- 11 • • • J- 



Jz(Ar,Aji,...,>)= / n^^' (E^^j 7 X N-D/2 (103) 

1=1 j=i (-T-S-T/2-i6 



The denominator becomes singular if the quadratic form Q{t) = —T ■ S ■T/2 approaches zero 
inside the integration region, explicitly: 



N~l N-1 

Q^^^ = 2^ ^''^'^^ + 5Z ^i^J- + Z = 0, (104) 



j,k=l j=l 
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where Xj^, Yj and Z are defined through Sjk- The singularity is protected by the i6 prescription 
or a finite width in some propagators. Viewing the integration volume as an (A^— l)-dimensional 
hypercontour in a space with — 1 complex dimensions, we are looking for a contour defor- 
mation which leads to a smooth and bounded integrand without intersecting the singularity 
hypersurface defined by eq. ()1()4|1 . For the analytic continuation of Q{i) — > Q{x) we make the 
ansatz jSE] x = t — i f. Now: 



Q(x) = Q{t)-^Yl ^JkTjTk - ^ 5Z ^ {X.ktj + Yk) (105) 

j,k=l k=l j=l 

This suggests the following choice for the deformation vector r: 
x{i) = t- i T{i) 

N-1 

^ ^ Ki^-tkr EiXjktj + Y,) ifr, = o ^^Qg^ 

' if Tfc ^ 



Note that the implementation of such a contour deformation also for the parameters which 
correspond to a non-vanishing width would not lead to the desired result, as lm{Q) < can 
not be guaranteed then. For X, a, (3 > the deformation moves the integration contour away 
from the poles, i.e. Im(Q) < 0, without causing any harm at the boundaries. 

The invariance under diffeomorphisms of the contour, Cx, means 

/ lldxif{x)= / Hdxifix) (107) 
1=1 J 1=1 



Co 



or in the given parametrisation: 

/ n = / n det ( fit- irit)) (108) 

'-^ '-^ 

The Jacobian is defined by 

-^ = 5i,-tX tf-\l - tif-^ [5i^ [a{l - ti) - (3ti] Xikh + Y^ +ti{l- ti)Xi^ 5{Ti) 



k=l 



Here 5{Ti) is equal to one if = and equal to zero else. To prove the invariance one can 
closely follow the derivation presented in the appendix of pni; apart from the fact that in our 
case a surface term is present which however turns out to be zero for the proposed contour 
deformation ()106|) . 

Some comments are in order. To simplify the discussion let us assume that F^ = for all I. 
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While AV ■ Q controls the size of the deformation, a, [3 control the smoothness of the de- 
formation at the integration boundaries. The vanishing of the gradient of Q{x) inside the 
integration volume is only critical if Q ^ at the same time. Using X^^ = X-^det{X), 
this critical situation is defined by the equations 



N-l 



det(X) Z=IY. ^u'^^^^ 



2 ^--y -'-^ 

N-l 



0<-sgn(det(X))5^X,7iyi<|det(X)| , j e {1, . . . , iV - 1} (109) 



1=1 



which are polynomial in the kinematical invariants. This situation corresponds to an 
anomalous threshold which is an exceptional kinematical configuration. With respect to 
integration over the phase space of external particles this is an integrable singularity. The 
critical surface in the integration regions is given by 



N-l 



t, = tf^*-=-$:x/F, (110) 



/=1 



All one has to do is to split the integration hypercube subsequently into 2^ ^ parts. This 
maps the critical surface to the integration boundary where adaptive integration routines 
can cope with the problem. If t^^^^- = or 1 subleading singularities may be probed. 

In the case of an UV divergent integral it is easy to explicitly isolate the UV pole. At 
one loop, only overall UV divergences are present, which manifest themselves in terms of 
F-functions in front of the parameter integral. With D = 4 + 2m — 2e one has an UV 
divergence, if and only if iV < 2 + m. Then r(A^ — D /2) has a single pole in e. In this 
case we need the sector integral Ji to order e : 



1 



N—^ N 

r 'L^ / " N N-4-2m-r , \ n 2+m-iV 

i 1=1 j=i 

1 - elog(--T . 5 ■ T - z5) + 2elog(5^T,) + ^(e^)] (111) 



X (T,,...T,v 



U N < 2+m, this integral is numerically stable without any modifications. The 0{1) term 
is a real number multiplying the UV pole. If X = 2 + m, there is a logarithmic singularity 
in the 0{e) term which can be dealt with by using the same contour deformation as 
defined in eq. (IKHij) . Note that no IR poles are present in the case N < 2 + m (the UV 
case). The UV and IR problems are thus nicely separated. 

The case of IR divergent 3-point functions can be treated analytically^. We give a com- 
plete list of 3-point functions for internal masses equal to zero in appendix IHl 



^It is also possible to extract the infrared poles by the method of iterated sector decomposition 82 from the 
parameter representations. This can be done in a completely automated way. A Laurent series in e is produced, 
where the coefficients are again Feynman parameter integrals with slightly more complicated denominators. 
The contour deformation can then be applied to a reduced number of Feynman parameters. 
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We note that the presented method should also be applicable for the direct numerical com- 
putation of finite Feynman parameter integrals without doing any algebraic tensor reduction. 
This is presently under investigation. 



7.2 Comparison between numerical and algebraic approach 

We present now a comparison between the algebraical and the numerical approach of evaluating 
the basis functions of our reduction formalism. 

As explained above we have algebraic representations of the higher dimensional 4-point func- 
tions where inverse Gram determinants are present, see sectional In a realistic application one 
expects that these representations will be numerically well-behaved in the bulk of the phase 
space under consideration. However, when approaching exceptional kinematical situations, 
compensations of large numbers will happen which will finally spoil a reliable evaluation. 

On the other hand, one expects that the method where integrals with Feynman parameters in 
the numerator are evaluated numerically from the start should not be sensitive to the presence 
of inverse Gram determinants, as there are no singular denominators in this representation. We 
have implemented the method for {n + 2)-dimensional box and ra-dimensional triangle functions 
for general kinematics and obtain good numerical behaviour using standard deterministic and 
Monte-Carlo methods. The two- and three-dimensional integral representations allow for a 
reliable direct evaluation of the required integrals. Comparing our numerical implementation 
to the algebraic one, we found that the algebraic implementation is much faster and accurate 
in the interior of the phase space, while the numerical one allows for the automated evaluation 
of the integrals near exceptional momentum configurations. This will be illustrated in the 
following by a simple example. 

Let us consider a 2 ^ 2 process pi + P2 — * Ps + P4 with pf = P2 = M'^ , where the momenta 
are parametrised as 



The Gram determinant is given by det(G) = 32M^ (1 + x^)"^ sin^ 6'. Exceptional configura- 
tions are the forward/backward scattering region, 6 = 0,7t and the case x = 0. 

In Fig. 17.21 we plot the real and imaginary parts of the functions -^4(2:4) with the 

parameter x varied from 1 to 0. For the plots we have set M = 1 and 6 = Tvr/SO. The 
output for the plot was obtained by a Fortran code working in double precision. Whereas the 
6-dimensional box function with numerator equal to one, shows a perfect agreement with 

the result from the algebraic representation for x values as low as x ~ 10~^, the box with a 
numerator, /|(^4), already starts to fluctuate severely for values of x as large as x ~ 10"^. 
However, for both of these integrals the numerical instabilities in the algebraic implementation 
occur in a region where the result could be safely extrapolated to the boundary a; = for 
our kinematical situation. In Fig. 17.21 the same plots are shown for the functions Iflz^z^) and 
/I (^3^1). The algebraic representations of these integrals have higher powers of inverse Gram 
determinants and are thus less stable. It is important to note that the instability occurs already 




Pi 



{E{x),0,0,Mx) , 
E{x) (l,0,sin^,cos^) 



P2 = {E{x),0,0,-Mx) 

Pi = E{x) (1,0, — sin 9, — cos 9) 



MVl + x^ 
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Figure 3: Real and imaginary parts of the basis integrals /4(1) and 1^(24), plotted versus 
the parameter x which interpolates between exceptional and non-exceptional kinematics, as 
explained in the text. The solid line stems from the numerical implementation, the dashed 
curves show the numerical behaviour of the algebraic representation. 

for values of x where one cannot yet safely extrapolate to the integration boundary. We have 
checked that applying quadruple precision for the evaluation of the discussed integrals leads to 
an improved behaviour. The function can be evaluated correctly for the whole plotted 

x-range for the given kinematics. The functions /|(^4), iKz^Zi) and If^z^zl) show instabilities 
below X ~ 10^'^, X ~ 10"^ and x ~ 10~^ respectively. The numerical problems are confined to 
a smaller phase space region but are still present. Of course the analytical approach could be 
improved by an algebraic expansion of the expression around critical regions, as has been done 
for example in 133] . Methods relying on such a Taylor expansion are faster, but require 
additional manual work, whereas the method suggested here is automated. 

Note that in all cases the purely numerical implementation is completely stable. It is actually 
possible to evaluate the given integrals numerically for all degenerate cases x = and 6 = 0,7t. 

The evaluation time for each plot point with an accuracy of better than one per cent using the 
numerical method is of the order of seconds on a standard PC with a Pentium 4 processor if 
Monte Carlo methods are applied. A precision of one per cent for the higher order correction 
should be well sufficient for phenomenological applications. The analytical evaluation is of the 
order of milliseconds and precise to standard Fortran double precision. Given the fact that the 
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Figure 4: Same as Fig. 17.21 but for the basis integrals iKz^z^) and Ifizszl). 

numerical evaluation is only called in a very small fraction of the phase space the relative speed 
is compensated by the small size of the critical regions. 

We conclude from this study that the best procedure for practical applications is to use nu- 
merical implementations near exceptional kinematical configurations and analytical ones in the 
interior phase space domains. 

8 Recipe for the practitioner 

In this section we would like to summarise how to apply the main results of the article in a 
practical calculation. Two steps have to be distinguished: first, expressing an amplitude in 
terms of our basis integrals and second, the evaluation of the basis integrals. 

Expressing the amplitude in terms of basis functions 

After having generated the amplitude as a combination of Feynman diagrams, it has to be 
expressed in terms of tensor integrals as defined in eq. (jl4j) . The further processing of these 
tensor integrals depends on the number N of external legs. 

In Fig. El we show the decision tree which indicates where to find the formulae to reduce an 
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A^-point tensor integral to our basis integrals. The latter are 2-point functions, n- and {n + 2)- 
dimensional 3-point functions and (n + 2)- and {n + 4)-dimensional 4-point functions with 
up to three Feynman parameters in the numerator. To perform the reduction with algebraic 
programs it is sufficient to look up and code the given equation numbers. Up to this point 
the representation is free from inverse Gram determinants. Algebraic simplifications might be 
applied after reduction for the coefficients of a given basis integral before proceeding to the 
numerical evaluation of the amplitude. We give a list of useful relations in appendices O and 

E 

As an illustration for a tensor reduction to basis integrals, let us show two examples. The first 
is the explicit expression for a rank two 5-point integral: 



dk — , 



It. E [Siis^'^-\H + s;,\si^}^\,,]ius\{j,k}). 



2 

jes kes\{j} 

As one can see, there is no inverse Gram determinant, only the inverse of the kinematic matrix 
S is present^. 

A rank one 6-point integral has the form 

I',''''ia;S) = -J2c^JUS\{j}) (112) 



j,ies kes 



m&S\{j,k} 



In both examples the basis integrals are already scalar integrals without Feynman parameters 
in the numerator. For higher rank tensor integrals this is not the case anymore. 



Evaluation of the basis integrals 

The case N = 2 needs no extra discussion. All necessary formulae are gathered in appendix 
1X1 For the evaluation of the remaining basis integrals we distinguish two cases, depending 

®The determinant detiS is vanishing only if one hits an anomalous threshold |80[ I5T] . which corresponds to 
an integrable singularity. The location of anomalous thresholds depends on the external kinematics and the 
internal particle masses. Note that there are no anomalous thresholds in the physical region if all internal 
masses are zero. This is also true for the massive case as long as the external particles are stable. 
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on whether the Gram determinants are numerically problematic or not. By eq. (j45|) Gram 
determinants are related to the quantity B, the sum of reduction coefficients. Technically, the 
splitting into safe and problematic regions can be achieved by introducing an adequate energy 
scale A. This scale should be chosen such that for BA^ > 1 the evaluation of the analytical 
expression for the basis integrals in terms of purely scalar integrals is numerically stable. In 
this case the evaluation of the basis integrals for N = 3 can be done by analytic reduction 
to scalar integrals by using eqs. (jH^ - lffOjl and eqs. (jA.ljl - fjA.lOj) . If the 3-point functions are 
IR divergent, reduction formulae are not applicable as detS = 0. In this case one has to use 
appendix El where explicit representations for all 3-point functions with and without Feynman 
parameters in the numerator for massless internal particles are provided. We do not quote the 
formulae for the IR divergent integrals where internal masses are present. Note that all 3-point 
functions with IR poles are relatively simple functions which can always be treated analytically. 

We illustrate our evaluation strategy for the 3-point functions in Fig. IHl where we indicate all 
relevant sections or equation numbers needed for the implementation of our method. 

The evaluation strategy for the = 4 basis integrals is depicted in Fig. [7| There are no IR 
divergences in this case. Note that the analytic branch in Fig. [3 contains implicitly evaluations 
of 3-point functions given in fig. IHl and an evaluation of scalar box integrals in 6-dimensions. 
For massless propagators analytical representations can be found in [IHl EI]- To our best 
knowledge no complete list of these integrals for all combinations of internal and external masses 
is available in the literature. However, we emphasise that the numerical evaluation, using the 
contour deformation method introduced in section 17. H can always be used if the analytical 
representation is not known. Analytical and numerical representation are complementary to 
each other in the sense that the former are fast and accurate in the main part of the phase space. 
The slower but robust numerical implementation also works well for exceptional kinematical 
configurations. As these dangerous phase space regions only cover a small part of the phase 
space, the speed issue does not pose a problem. 
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Figure 5: Reduction of A^-point tensor integrals to basis integrals. The indicated equations can 
be implemented directly into an algebraic computer program. /]v(l|ii|ii;i2|ji, j2, js) denotes 
the integral with zero, one, two or three Feynman parameters in the numerator. 



35 




Figure 6: Evaluation of the basis integrals: the triangle case, 
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Figure 7: Evaluation of the basis integrals: the box case, N = 4. For the analytical evaluation 
of eqs. use Fig. IHl 
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9 Conclusion 



In this paper, wc have presented a complete method for the calculation of one-loop nmlti- 
leg amplitudes. In principle it can be applied to arbitrary A^-point problems, the limitation 
coming only from computer power. We offer a new method for tensor reduction which has 
several advantages: 

• The formalism, using dimensional regularisation, is valid for both massless and massive 
particles. 

• Infrared divergences are easily isolated by construction. They appear only in terms of 
three-point integrals which we list explicitly for the massless case. 

• The formalism is completely shift invariant, such that its iterated application does not 
require the redefinition of loop momenta. 

• Integrals in more than n = 4 — 2e dimensions do not have to be evaluated for > 5 
external legs, as we have proven that they drop out. 

• By using Feynman parameter integrals with non-trivial numerators as basis functions, 
inverse Gram determinants can be completely avoided. 

We present two methods to compute the basis integrals which are the endpoints of our reduction. 
First we discuss the possibility of a "purely algebraic" approach, where all non-scalar integrals 
are reduced further to end up with scalar integrals only. This procedure re-introduces inverse 
Gram determinants, which could spoil the subsequent numerical evaluation if an exceptional 
kinematic configuration is approached. We also propose a method to compute the basis integrals 
completely numerically by multi-dimensional contour deformation in Feynman parameter space. 

We have compared the two alternatives, with special emphasis on the behaviour for exceptional 
kinematics. We show that the semi-numerical approach, where non-scalar integrals are used 
as basis integrals, is very stable if an exceptional kinematic configuration is approached. On 
the other hand, the evaluation of the same integrals in a representation where they have been 
reduced algebraically down to scalar integrals, is stable only in the interior of the phase space. 
In this region however - which is the bulk of the phase space - their evaluation is of course 
faster than evaluating the non-scalar form. In a program to calculate one-loop amplitudes, it 
is possible to combine the virtues of both alternatives by using the algebraic representation in 
the interior phase space domains and switching to the semi-numerical one at the phase space 
boundaries. 

The paper contains a complete list of form factors for integrals of rank r < N, where N = 
1, . . . , 5, and it is shown that for iV > 5 it is not necessary to introduce new form factors. 
The formalism naturally maps tensor A'^-point integrals with N >6 to combinations of 5-point 
integrals and reduction coefficients. 

All the formulae needed for a direct implementation of the formalism are given in the paper, 
except some well-known 3-point and 4-point integrals, such that readers who are less interested 
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in the technical details can straightforwardly use the method by following the guidelines in the 
section "recipe for the practitioner" . 

Further, we list several relations, between form factors as well as between reduction coefficients, 
which can be very useful to achieve a more compact form of a given amplitude and/or to perform 
checks of the program. In particular, we give a helicity decomposition for the case of massless 
6-point functions which leads to very compact expressions for the reduction coefficients. 

In summary, the formalism presented here, as it can deal with massive particles as well as 
infrared divergences, has no restriction on the number of external legs and is suitable for 
numerical integration, can be used as the basis of a general program to calculate multi-leg 
one-loop amplitudes efficiently in a highly automated way. 
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A Form factors for A^^ = 1, 2 

For completeness, wc provide the one- and two-point functions in this appendix. We give the 
kinematical arguments here in terms of invariants, as S is trivial in these cases. 

ir'(a,;ml) = AZj^mj) ^ A^^mj 1^(0,0, mj) (A.l) 

Here we use the fact that 1-point functions may be written as degenerate 2-point functions. 
The kinematical matrix for the 2-point functions is 

\ —S + m( + 1712 2 777,2 / 

The Lorentz tensor decomposition for the tensor 2-point functions 

I-''^\ar;s,mlml) = J] Af^^^ 4'^(5) (A.3) 

hes 



39 



defines the form factors 



{S, 171^,1712) 



B'^\s,mlmi) 



1 



m: 



2 mlll^{s, ml ml) + m^/^^lO, 0, ml) 



2s 



2{n - 1) 



-S + — 777,2 



I^{s,mlml) 



7777 



777 



Aif(s,mi, 7772) 



2(77 - l)s 



2 [/2"(s,7n?,m^)-/2"(0, 777^, 777^)] 
/2 (5,7771,7772) 



77(5 — 777^ + ml] 



m^ 



m 



'mis^mlmi) - r^{0,mlmi)] 



-2 I^{s, ml ml) + (77 - 2) 777f/2''(0, 0, m{ 
The other tensor coefficients are obtained by the following relations 



Al^ + AY 

42,2 I a2,2 



12 
2,2 
22 

^2,2 I ^2,2 I <2,2 I ^2,2 



HI 
i2.2 
^21 
.2.2 
^21 



Al?- + A: 



-'2 



(A.4) 



(A.5) 



(A.6) 



which follow directly from the fact that in one-loop parameter integrals, the sum of all Feynman 
parameters is equal to one. The general 2-point scalar integral is well known [7S] ■ We give here 
the integral representation for completeness. 

1 

12(3, mi, m2) = T{e)- j dx\og{-sx{l-x)+xml + {l-x)ml-iS) + 0{e) (A. 7) 



If the external vector is light-like the formulae degenerate to 

(4 mf — n m\ + 77 ml) m\ 



A{ (0,mi,'7772 



52'2(0, 777^,777^ 



+ 



2r7(r772 



777 



2\2 



/2 (0,0,777f) 



^ 777^^ - nmf-F 77 777^)7771 ^^^^^^ 



2n{m\ 



777f)2 



^ 2T K ^2^(0, 0, ml) - ml 1^(0, 0, 777^)) 



777T 



77(7773 ^ ^1) 



1^(0,0, ml + 



I-iO,0,ml) 



n{m\ — ml) 
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-Ami {ml -ml) -2 ml)] ml 12(0,0, ml) 

-8 ml 1^(0,0, ml)) (A.8) 



where 



l^{0,ml,ml] 
l^{0,0,ml] 



mllUO,0,ml)-mll^iO,0,ml] 



Il'{ml) 



2 2 
m2 — m^ 



n 



-mO,ml,ml) = -T{l-n/2) Kf/^- 



(A.9) 



In the case mi = m2 one finds 



Al'\0,ml,ml) = !^I-(0,0,ml) 
B''\0,ml,ml) = ^/,"(0,0,m?) 

Al',\0,ml,ml) = !^I-(0,0,ml) (A.IO) 



Two-point functions with no scale at all are defined as zero in dimensional regularisation. 



B Divergent three— point functions with massless prop- 
agators 

We consider here only three-point functions with massless propagators. The parameter repre- 
sentation is given by 

mJuJ2,J3;S) = -T{3-^) ffldz,6{l-J2^i)^n^n^^siR')^-' (B-H) 

i=l 1=1 

= -Z1Z2S12- Z2Z3S23- ZiZ-iSis -i5 

If one or two invariants out of the set {Su, S23, vanish, one gets an IR divergence. In this 
case one has det5 = 0, such that the formulae given in section I^TTl do not apply. Therefore we 
provide analytic representations for all 3-point parameter integrals with massless propagators 
here. An overall coefficient r-p is defined as 

r(i + e)r^(i-e) 

^ r(l-2e) 

For three-point functions with one non-zero invariant, denoted by X, we labelled the inter- 
nal propagators in such way that 1S13 = X and 1S12 = S23 = 0. For three-point functions 
with two non-zero invariants X and Y, we set 1S23 = X and S13 = Y. Thus the integrals 
I^{zi, Zj, . . . ; 0, 0, X) are symmetric under exchange of Zi ^ Z3 and the I^^Zi, Zj, . . . -,0, X, Y) 
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are symmetric under simultaneous exchange of Zi Z2 and X ^ Y. We obtain for integrals 
with no Feynman parameters in the numerator: 



IU0,X,Y) 

with one Feynman parameter: 

mi;Q,Q,X) = 
/3-(2;0,0,X) = 
I^{3;0,0,X) = 

I^{1;0,X,Y) = 
I^{2;0,X,Y) = 
I^{3;0,X,Y) = 
with two Feynman parameters: 

/3-(l,l;0,0,X 

/3-(2,2;0,0,X 

I^{3,3;0,0,X 

I^{1,2;0,0,X 

/3"(1,3;0,0,X 

I^{2,3;0,0,X 



"^HoiX, -e) 
!lH,{X,Y,-e) 



rr 



€ l-2e 

rr 1 

^ l-2e 
rr 1 



e l-2e 

rr 1 — e 

^ l-2e 
rr 1 — e 



Ho{X, -e) 



rr 1 



e l-2e 



^2(^,1^, -e) 
H2{Y,X,-e) 
H,{X,Y,-e) 



I^{1,1;0,X,Y 
I^{2,2;0,X,Y 
IS{3,3;0,X,Y 
I^{1,2;Q,X,Y 



rr 1 
e 2(1 -2e) 
rr 1 



e2 (l-e)(l-2e) 
rr 1 



Ho{X, -e) 



e 2(1 -2e) 
rr 1 



e 2(l-e)(l-2e) 

rr 

Ho{X, -e 



2(l-6)(l-2e) 
rr 1 



e 2(l-e)(l-2e) 
rr 2 - e 



HoiX, -e) 



e2 2(1 - 2e) 
rr 2 - e 
^ 2(1 - 2e) 

rr 1 



^3(^,^,-e) 
H,{Y,X, -e) 
i7i(X,F,-6) 



e 2(1 -2e) 



B.12) 
B.13) 

B.14) 
B.15) 
B.16) 

B.17) 
B.18) 
B.19) 

B.20) 

B.21) 

B.22) 

B.23) 
B.24) 

B.25) 

B.26) 
B.27) 
B.28) 
B.29) 
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/-(l,3;0,X,r) = 
/3"(2,3;0,X,F) = 

with three Feynman parameters: 



" e 2(1 - 2e) 
rr 1 
" e 2(1 - 2e) 



H2{X,Y,-e) 
H2iY,X, -e) 



1^(2,2,2 
/a (3, 3, 3 

73^(1,2,2 
/3(1,1,3 
/3(2,2,3 
^"(1,3,3 
/3(2,3,3 
/3(1,2,3 



0,0, X 
0,0, X 
0,0, X 
0,0, X 
0,0, X 
0,0, X 
0,0, X 
0,0, X 
0,0, X 
0,0, X 



/3"(1,1,1;0,X,F 
/3"(2,2,2;0,X,F 
73"(3,3,3;0,X,r 
73^(1,1, 2; 0,X,r 
73^(1, 2,2; 0,X, r 
73^(1, 1,3; 0,X, r 
73"(2,2,3;0,X,r 



rr 



2-e 



e 2(3-2e)(l-2e) 
rr 3 



e2 (l-e)(3-2e)(l-2e) 
rr 2-e 



Ho{X, -e) 
Ho{X, -e) 



e 2(3-2e)(l-2e) 
rr 1 
T2(3-2e)(l-2e) 
rr 1 



e (l-e)(3-2e)(l-2e) 
... . . ^ HoiX, -e) 



77o(X,-e) 
77o(X,-e) 

77o(X, -e) 



2(3-2e)(l-2e) 
rr 1 



e (1 -e)(3-2e)(l -2e) 
... . ^ Ho{X, -e) 



Ho{X, -e) 



2(3-2e)(l-2e) 
rr 1 



Ho{X,-e) 



€ 2(3-2e)(l-2e) 

^^ 2(3-2e)(l-2e) ''°^'''-^^ 
rr (3-e)(2-e) 



e2 2(3-2e)(l -2e) 

rr (3-6)(2-6) 
£2 2(3-2e)(l -2e) 
rr (2 - e) 



e 2(3-2e)(l-2e) 
rr (3-e)(2-e) 



774(X,F,-e) 
774(y,X,-e) 
H,{X,Y,-e) 



e2 2(3-2e)(l -2e) 
rr (3-e)(2-e) 
e2 2(3-2e)(l -2e) 
rr (2 - r) 
e 2(3-2e)(l-2e) 
rr (2 - 6) 
e 2(3-2e)(l-2e) 



(773(X,y,-e)- 774(X,y,-e)) 
(773(y,X,-e)- 774(r,X,-6)) 
HsiX,Y,-e) 
HsiY,X,-e) 



(B.30) 

(B.31) 
(B.32) 



(B.33) 
(B.34) 
(B.35) 
(B.36) 
(B.37) 
(B.38) 
(B.39) 
(B.40) 
(B.41) 
(B.42) 

(B.43) 
(B.44) 
(B.45) 
(B.46) 
(B.47) 
(B.48) 
(B.49) 
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I^{1,3,3;0,X,Y) = ""^ 

I^{2,3,3;Q,X,Y) = 

I^{1,2,3;Q,X,Y) = 
Higher dimensional 3-point integrals: 

IS+\0,X,Y)^ 



(1 



e 


2(3 


-26)(1- 


26) 


rr 




(1-^) 




e 


2(3 


-26)(1- 


26) 


rr 




(2-6) 




e 


2(3 


-2e)(l- 


26) 



H2{X,Y,-e) 
H2{Y,X,-e) 

{H2iX,Y,-e)-Hs{X,Y,-e)) 



rr 



1 



6 2(l-6)(l-26) 

rr 1 



i^o(X, 1 - 6 



73"+2(l;0,0,X) = 



/3"+^(2;0,0,X) = 



rr 



6 2(l-6)(l-26) 

rr 1 
T2(3-26)(l -26) 
1 



/3"+2(3;0,0,X) 

/3-+2(l;0,X,F) = ''^ 
/3"+2(2;0,X,F) = 



6 2(1 - 6) (3 - 26) (1 - 26) 

rr 1 



6 2(3 - 26) (1 - 26) 

(2-6) 



H,{X,Y,l-e 

Ho{X, 1 - 6 
i^o(^, 1 - e 
i^o(^,l-e 



6 2(l-6)(3-26)(l-26) 

rr (2 - 6) 

6 2(l-6)(3-26)(l-26) 

rr 1 



7o"+2(3;0,X,y) 

^ ^ ' ' ' ^ 6 2(3 - 26)(l-26) 
The functions Hq, Hi, H2, and H^^ are given by: 

~Y 



H2{X,Y,l-e 
H2{Y,X,l-e 
H,iX,Y,l-e 



Ho{X,a) 
H,{X,Y,a) 
H2{X,Y,a) 
H3{X,Y,a) 
H^iX, Y, a) 



X-Y 



+ 



Y-X l + a {Y-Xf 



ya 2 

+ 



y2+Q; _ J^2+q; 



y-X 1 + Q;(y-X)2 (1 + q;)(2 + q;) (F - X)^ 

"p« 2 yl+« g y2+a 

r - X ^ 1 + a (r - X)2 ^ (1 + a) (2 + a) (F - X)3 

+ 



(1 + a) (2 + a) (3 + a) {Y - Xf 
X-i5 



(B.65) 



All these functions have a regular behaviour when X — Y. It is for this reason that we want 
to keep them as they are. If the coefficient in front of them is proportional to {X — Y), then 
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they can be reduced using the following properties: 



Hi{Y,X, 

H3iY,X, 
H,{Y,X, 





- X) Hi{X,Y,a) 


= Y — X 




(B.DDj 


(Y 


- X) H',(X,Y,a) 


1 + a 


^ X Hi (X, y, a) 

I + a 




(Y 


-X)H,{X,Y,a) 




2 

- XH2{X,Y,a) 
L -\- OL 


(B.68) 


(y 


- X) Hi{X,Y,a) 


_ ^ ya 

3 + a 


--J—XHs{X,Y,a) 
6 + a 


(B.69) 


a) 


= H,{X,Y,a) 








a) 


= H,{X,Y,a)- 


-H2{X,Y,a) 






a) 


= H3iX,Y,a)- 


-2H2iX,Y,a) 


+ Hi{X,Y,a) 




a) 


= -Hi{X,Y,a) 


+ 3 HsiX, Y,a)-3 H^iX, F, a) + H^{X, F, a) 


(B.70) 



We use n = 4 — 2e where e < in the infrared region. 

C Proof of the absence of higher dimensional integrals 
and Gram determinants for = 5 

In the next subsection, we will show that the coefficient multiplying the higher dimensional five- 
point integrals is of order e. We will repeatedly use the fact that for general 5-point kinematics, 
the metric tensor in 4 dimensions is expressible by a tensor product of external vectors. Many 
simplifications will occur by neglecting terms of "(9(e)". For scalar quantities it is clear what 
that means. For tensors, we say that a tensorial structure is of 0{e) if differences of tensors 
defined in n and in 4 dimensions occur. Contracting such differences with kinematical objects 
like external momenta, polarisation vectors or fermion currents will always lead finally to scalar 
quantities of 0{e), which can be neglected in phenomenological applications at one loop. 

The price to pay for the disappearance of the higher dimensional integrals is that inverse 
Gram determinants (1/-B) are reintroduced explicitly. In subsection IC.21 we show how these 
spurious divergences cancel out analytically. This will lead us to a representation of 5-point 
functions which is free from higher dimensional 5-point integrals and 1/B terms. The corre- 
sponding form factors are listed in section IHl of the main text. 

C.l The fate of higher dimensional integrals 

For N = 5, the following relation, shown in appendix ID.3| 

= ^ 

will enable us to remove the higher dimensional five-point integrals 
are defined in eqs. (|HT|l and 
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(C.71) 

. The tensors and T^^ 



For rank zero and rank one, it is trivial to see that the coefficient of I^'^'^ is of order e. For rank 
0, one has, according to eqs. and 



I-{S) =Y.hI2{S \ {j}) - (4 - n)BI-^\S) . (C.72) 

For rank 1, one obtains by apphcation of ()46|) : 

^5 '"(a, 5) = - E ^4(5 \ {j}) + (4 - n) ir\S) . (C.73) 

For higher rank (r > 2), we prove by induction on r that the coefficient of the higher dimensional 
five-point integrals is of order e. 

In step one we show that the assumption is true for r = 2. Indeed, in the case of rank 2, direct 
application of (jlUI) leads to 



I^'^\ai,a2]S) = I dk ^^^^^^ 



-E^.^a.^4'>i;^\{j})- (C.74) 

The integral ; S) with Feynman parameter Zi in the numerator can be expressed in terms 

of I^~^'^{S). For this purpose, we contract both sides of eq. ()C.74|) with A'^^, where b and c are 
arbitrary labels in S. Using 

qH^bc, = l{ql-ml-[ql-ml]-Sa, + Sac) (C.75) 
r,7A,rf^ = {Sac- Sad) (C.76) 

V(^A,rf^ = ^ (Sac -Sad) (C.77) 

C^a^cdf. = I {Sjd-S,c + b,[Sac-Sad]) (C.78) 

leads to 

j2ir\^;s)A^^, = ^\-v^air\s) 

E ■ (C.79) 

The insertion of the latter into eq. ()(I74j) yields 

I-''"'{a„a,;S) = S^-lr^a, + {n - 3)^^^^ It\S) 
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+ ^ E ^4 \ {^}) • (C.80) 



B 



Using now eq. ()(173|1 for the last term in eq. ()(I80j) results in 

(C.81) 

Now we can use eq. fjC71|) to see that the coefficient of in eq. ()C.81|) is indeed of order e. 

Let us now assume that for rank r — 1, all the higher dimensional integrals involve an 0{e) 
tensor of the type T'^^^ — 2 V(^^ V^^Z-B, and show that this is also the case for rank r. Firstly, 
we carry out the momentum integration in the first term of eq. (jlUj) . using the shift (j^^ and 
eqs. (jm, (jS2I)- This leads to^: 



+ dZ dl 



N 



where 



N 

dZ = T{N) l[dzJ{l-Y,^i) 

ieS 1=1 

and the qa denote the ga- vectors in terms of the shifted loop momentum given by eq. ()47|) . 

Secondly, we contract eq. ()C.82|) with A^c^^, where b,c E S are arbitrary, using ()C75|) for 
the left-hand side and ()C.76|) to ()C.78|) for the right-hand side. This yields 

C • • • C'-"i V- , f rf ((ij - €1 ■ ■ ■ 



Using now eq. ()(y.83|l to replace the term containing J2ies 

A^. in jniHl, we obtain: 



jes ^ ^ 



^Here we keep N arbitrary to show that the relations derived are not only valid for N — 5. 
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' ar Tn,fj,i...^, — 1 



B 

+ / dZ 



dl 



ai, . . . , a.r_i ; S) 

I, (t;:;,^ - 2 vzvyB) . . . g^- 



AT 



. (C.84) 



Specifying = 5, we can see from eq. ()C84|) that a rank r 5-point integral can be written 
as the sum of 4-point integrals of rank r — 1 plus a 5-point integral of rank r — 1 plus a 
term which potentially generates higher dimensional 5-point integrals but is proportional to 

A few comments are in order here. 



(i) Eq. ()C.83|1 defines another way of reducing the tensor integrals, by which a rank r Ap- 
point integral is expressed as an infrared finite part plus a sum of rank r (A^ — 1) -point 
integrals. Contrarily to eq. ()C.82|) . eq. ()C.83|) reduces the number of propagators but not 
the rank. For this reason the way used in section |3] is preferable. 



(ii) In the case where S is not invertible, the derivation is equally valid, as eq. ()C.78|) still 
holds in this case. This can be easily seen by contracting eq. ()C.82|) for r = 1 with A^c^. 
One obtains 

c'.a^^c, r^-iis \ {j}) = \ {r^^As \ {c}) - r^^As \ m + {Sa, - r^{s)) . 

jes 

Now we can use I'^{S) = J2j<^s ^N-ii^\{j})^ which is always vahd if i? = J2jes ~ 0' 
and fulfilled to order e for N = 5. The reduction coefficients bj can be constructed as in 
section ITT| if S is not invertible. Thus one obtains 

J2 C';Ai>c, In-i{S \ {j}) = IY. ^N-iiS \ {j}) [S,c - S,, + h, {Sa, - 5..)] . (C.85) 

which leads to eq. ()C.78|) without using the inverse of S. 

(iii) A similar proof is given in eqs. (34) to (37) of jSHI; but now we have a shift invariant 
formulation, and all coefficients are expressed in terms of quantities containing only S~j^ 
instead of Hij (the pseudo-inverse of the Gram matrix G). 



In summary, we have shown that 



ig [tti, a2, ■ ■ ■ , ttr , O ) 



' ar 



^^^^ . . . , a,_i ; 5) + 5^ 6, /^'^-'^^-(ai. 



,ar-x\S\{jY) 



je5 



n,ii\--iir~\ ( ^ 

4 I'^l! 



ie5 

.ar-i;5'\{j}) 



(C.86) 
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In eq. ()C.86|) . inverse Gram determinants have been explicitly reintroduced in the term 



_^„,,....,„_. ^^^^ ^ ■,S) + J2 IT''-'-' (ai, . . . , ar.^ ; S \ {j}) 



(C.87) 



We show in the next section how they drop out. 



C.2 Cancellation of 1/B terms for = 5 



In this section we will prove by induction that F^^^^^^''{al, . . . ,ar), defined in eq. ()(187j) . is in 
fact free from 1/B terms if the factor V'^^/B is combined with the expressions inside the square 
bracket. 



To this end, we first use eq. ()C.86|) for rank r — 1 



1---Mr-1 



[ai, . . . , a._i) - CtaZ. (ai, • • • , «.-2 ; S \ {j}) (C.88) 

ie5 



and insert the above equation into eq. ()C.87|) . leading to 



~B~ 



0-1, ■ ■ ■ , Ctr-l ; 



4 [ai, ■ ■ ■ , 



j€S 



a._2 ; 5 \ {j}) + 5^ 6, rr'-'-\a,, . . . , a._i ; 5 \ {j}) 
flr-i ; S \ {]]) 



(C.89) 



Actually, the iteration of eq. ()C.86|) performed in eq. ()U.89|) singles out the pair of indices 
(a^-i, /ir-i), which hides the manifest symmetry of 1^'' '"^"^ {cli, ■ ■ ■ , ctr-i ; "S*) with respect to 
all pairs (ai,/ii), ■ ■ ■, (0^—1, /ir-i)- However, the exphcit cancellation of 1/B terms relies on 
this symmetry. Therefore we introduce a symmetrisation operator which acts as follows: If 
y^^ii...^is . . . , tts) is a tensor which already is symmetric with respect to the s — 1 first indices, 
'Eg is defined by 



S,[iy^-'^»(ai,...,a,)] = -[iy^-'^^(ai,...,a,)+c. p. 



(C.90) 



where "c. p." means the sum over cyclic permutations of (ai,/ii^ 
write eq. ()C.89|) as 



(as,/is)- Thus we can 



ig {ai, a2, ■ ■ ■ , ttr , o ) 



r-iF'^i-'^'-nai,---,a,_i) + Q 



fii...fj., — 1 



- E ^Tar • • • , ; 5 \ {j}) , (C.91) 

i65 
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where 



3&S 



C;;;!,/r'"-^(ai,...,a.-2;^\{j}) 



+ 6,./r^-^-^(ai,...,a._i;5\{j}) 



(C.92) 



In order to show that the 1/ B terms in the first hne of eq. ()C.91|) drop out, let us first consider 
Q^,^...^lr-l ^ Using eq. (in83|) with = 4 for the first term and eq. (|n82|l with = 4 for the 
second term in eq. ()(I92j) . we obtain: 



(C.93) 



(g| - m|) {ql - ml) q^l . . . q^;-_ 



where 



Ml-'-A'r-l 



h Kzi: - 2 c%-\ v:^ ) jdz jdi "Jfr (c.94) 



{P - R] 



2^4 



Q 



fJ,l...fJ,r-2 

2j 



Ri 



{P - R]Y 



(C.95) 



i=l 



k,ies\{j} 



and the identities (ID.145jl . (ID.146jl and (ID.147jl have been used. The last term of eq. (1(193^ 
vanishes due to antisymmetry with respect to j and k. Therefore, the first line in eq. ()C.9ip 
can be rewritten as 



' Or 



■^J■r-2-\ 



(C.96) 



The term multiplying Q^J-'''-' involves {bj -BCf^^) V^l/B, and due to eq. ^(TTH we have 

(C.97) 



a, 



explicitly free from 1/B terms. 



For the remaining contribution, we will show that Q^^'"^'' ^ = F^'^' '^'— ^ + 0{e). More in 
detail, we will prove by induction that 



pfll...^ir-l — \ ^ "'•-1 



i65 



r-2 
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and show by direct calculation that we also have 



/il.../tr-l 



-r-1 



V 



Mr-1 



r-2 



J ar-2 / "i-^J 



0(e) (C.99) 



Eq. ()(I99j) is established in subsection 1(131 for r = 1, . . . , 5. The case r > 5 will never be 
needed, but the mechanism is exactly the same. 



Let us now show eq. ()C98|) by induction. The induction start is r = 2, because for r = 1 the 
absence of 1/5 terms is trivial. For r = 2, we combine eqs. ()C.73|) and ()C.88|) to obtain 



F^^{a,) = {A-n)V^,lI-+\S) = 0{e) 



From eq. ()C.94|) . we obtain 



(C.IOO) 



(C.lOl) 



Now let us assume that ()(I98j) is fulfilled for rank r — 1. To prove the step r — 1 ^ r, we 
use eq. (1(19111 for 

; S) and replace F^^^-^''-^ by ()(".98p . which is true by the 

induction assumption, to obtain 



I^''^^'"^'' (cbi, ■ ■ ■ ,ar; S) 



a, 



f (- |E ^ ^r-. [{h VZzl - BC^-l) Q--— ] I 



'2i 



(C.102) 



Comparing with eq. ()(I99jl . we see that the term Sf—i {. . .} in the first line of eq. ()(I102j) 
is equal to Qi^ "^'^"^ Comparing the remaining terms to eqs. ()C.88|) and ()C.98|) proves our 
assumption. Therefore, using eq. ()C.97|) . we see that rank r 5-point integrals can be written as 



Tn,^I.l...^lr ( 



- C^a. Il'^'-^-' (ai, • • • , a._i ; S \ {j}) 1 + 0{e) 



(C.103) 



which is a combination of 4-point integrals. 



C.3 Auxiliary relations 

We do not have to compute 5-point integrals with rank bigger than five. Therefore it is sufficient 
to establish eq. ()C.99|) by direct calculation of QMi-^r-i g^j^^ Q'^^'^'"^ for r = 1, . . . , 5. Using 
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the definition (ITTQill and eqs. (ID. 1451) . (ID.14(ijl and (ID.147jl . one gets: 

(C.105) 



Q 



Ml M2 MS ^4 
1 



4 ^ 

I „A»2At3 r)/il/i4 _|_ „/i2A'4 r)A»lM3 _|_ ^M3M4 r)Ml/i2\ 
y ai a4 ' y ai as ' y a\ 02 J 

- J2 ir\hk;S\{j}) 

i,k£S 

X f A/^l A/^Z r)At3 At4 _|_ A Ml AM3 r)M2M4_|_AMl AM4 ^)^l2^l3 

\ ai i 02 fc 113 04 ai « as A: a2 04 ^ ai i 04 A; 02 03 

1/^/^2 AM3 M4 _|_ A M2 AM4 r)Ml M3 _|_ A M3 AM4 r)MlM2N 

~ a2 i 03 A: ai 04 02 i 04 cii 013 ' 03 j a4 A ai 02 / 



IM2 M3 

12 13 



where 



nMi M2 

ai 02 



aia2 '^ja2 ^ ai j ai "^a; 



2^ 
02 > 



Similarly, eq. ()C.95|) leads to 



Q% = {2-n)Y,K\^n''\^■.S\{J}) 



i&S 



Q 



Ml M2 
2j 



iG5 



qMi M2 M3 



'2i 



- E Aa,\A::2\A;^jr^(z,fc,/;s\{j}) 



(C.106) 

(C.107) 
(C.108) 



(C.109) 



(C.llO) 



If we insert now eqs. fjC.108|) to ()C110|) for Q2j into the right-hand side of eq. ()C.99|) and use 
eq. (fn97jl . we obtain the expressions ()(;.1()4|1 to for Qi. 



D Useful relations 

In this appendix we give a collection of formulae which are useful if one wants to perform 
algebraic simplifications of loop amplitudes. The relations in subsection ID. II can also be very 
useful to perform checks on the implementation of the form factors in a computer program. 
Some of the relations are already given in the main text, but the purpose here is to list them 
for quick reference. 
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D.l Relations between the form factors 

The following identities have been used extensively to obtain the relations given below: 

1 



qa-Abc = ^{q^ -ml-[ql-ml]-Sab + Sac) 



(D.lll) 
(D.112) 



Using relations (jD.lllj) and ()D.112j) and multiplying with a vector A both the definition and 
the expression in terms of form factors of an integral, one finds relations between the form 
factors. 

D.1.1 Four-point functions 

For example, we have 



dk 



qa-A 



be 



1 



dk 



^ ^t'^{S) {Sic — <Slb + Sab — Sac) 



(D.113) 



Eq. (ID.113jl implies 



J2Af\S)iSu-Sib) = A'^\S\{b})-A''\S\{c}). 



(D.114) 
(D.115) 



One can proceed in the same way for all the four-point functions and finds for rank 2: 



hes 

{Si,c-Si,b) 



hes 



where the definition 



-At;\S) (D.116) 
+ b A^;\S \ {b}) - 4 , Al\S \ {c}) , (D. 117) 



1 - Sji 



Oifj = / 



(D.118) 



53 



has been used. For rank 3, one obtains: 

J2Bt'\S) = -B^'\S) (D.119) 



i&s 

J2Bf'\S){Su-Si,) = B'^'{S\{b})-B'^'{S\{c}) (D.120) 
les 

E^hlJ^) = -^tliS) (D.121) 
hes 

E ^Xis(S) (^isc - Si,,) = 2 {Si,, - Si,,) BffiS) + 2 {Si,, - Si,,) BffiS) (D.122) 
hes 

+ Si, b ^l, b ^fiij ("^ \ {^}) ~^llcSl2C iS\{c}) 

For rank 4: 

J2Bf;i{S) = -BffiS) (D.123) 

J2Bt;i{S){Si,,-Si,b) = 2(5,,e-5,,,)C^'^(5) 

+ Kf, BffiS \ {b}) - Si,c BffiS \ {c}) (D.124) 



J2^tUu(S) - -Atfi,i,iS) (D.125) 

/4GS 

(4G5 i=h,l2,l3 

+ Si,,,Si.^bSi,bAl\i,^iS\{b}) 

— Si, c Sl2 C S13 c 

D.1.2 Five-point functions 

Below we will use the definition 

= 2 - S.A (D.127) 



Rank 1: 



J24'\S) = (D.128) 
5^ = A''%S\{b})-A'''iS\{c}) (D.129) 



Rank 2: 



ie5 



J2A'ifi,iS)iSi,,-Si,b) = 2(5,,,-5,,,)5^'2(5) 

+ m) - KcAtfiS\{c}) (D.131) 
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Rank 3: 



J2 Bf'\S) = -B^'^S) (D.132) 

E4'M3(^) = -4i(^) (D-133) 

[B'^'{S\{b})-B'^'{S\{c})]+2 ^ (5.c-5.6)4f,,}\«(^) 

i=h,h 

{S\{c}) (D.134) 



Rank 4: 



E<'m3u(^) = -<'m3(^) (D-136) 

E ^3 u (s) + Hi^i. Bli {s) + m,,, sji {s) + m,,, Bf;l {s)] {s^^- ,) = 

Hi,i, [2 (5,3, - 5,,,) C'^\S) + \ {^}) - 4ci?f3"'(^ \ {c})] 

+i/w3 [2 (5,,, - 5,,,) C'^\S) + 6i,,Bff{S \ {b}) - 4c<'(5 \ {c})] 
+i^,,,3 [2 (5,,, - 5,,,) + 6i,,Bff{S \ {b}) - KcBtl\S \ {c})] 

+2 E ('^--'^^'>)<..u\»(^) 



+ 6 b Sis b ^tfh h ('^ \ {^}) - C 4 C 4 C ^t'l /3 ('^ \ i^)) 

(D.137) 



Rank 5: 



E<'m3uJ^) = -4'm3J^) (D-138) 
(-f^h h His k + -f^/i ^3 Bi2 u + -f^h U Hl2 la) E 

+Hh h E U h + ^'i 'a E ^hk h ('^) + k E ^?2l3 is ('^) 

ises Jbgs ises 

+/7,,,3 e^?;m5(^)+^'^'^ E^M3J^)+^'3u E<iJ^) 

iseS iseS h&S 

-{Hi^ Za -^^3 «4 + Hl^ (3 -ffzj (4 + Hi^ ^ if;2 ,3) C^''^{S) 

-Hh I2 -Sfjt (-S") - Hi^ Is B^^]^ {S) - Hi^ i^ Sjig (5") 
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5,4 



,5,4 



(D.139) 



+ Hh k 

+ Hh h 



2 (5,,, - 5,,,) (5) + 2 (5,,, - 5,,,) CffiS) + J2 i^i.c- Si,,) J5) 



2 (5,,, - 6i,,) ClfiS) + 2 (5,3, - 5,3,) Cf/(^) + 5^ {Si,,- Si,,) Bf;l,^{S) 



2 {61,,- 61,,) CffiS) + 2 {6u, - 5u,) C.f (S) + J2 (Skc- Si,,) Bf;l,^{S) 



2 {Si, , - 5,, e) Cf/(^) + 2 (5,3 , - 5,3 ,) CffiS) + J2 (Sk c - Si, ,) B_ 



+ Hi 



h h 



2(5ufe-5,,e)Cf (5) + 2(5, 



'}hb - 5i^c) Cf;\S) + {Si,c- Si,b) B 



+ J2iShc-Si,,)A',; 



h h h I. 



SS) 



+ 2 {5i, , - 5i, ,) BfX J5) + 2 (5,, , - 5,, ,) Bf;l (S) 
+ 2 (5,3 , - 5,3 ,) Bf;l J5) + 2 (5,, , - 5,, ,) Bf^ ,3 (S) 
{Hi, Hi, I, + Hi, ,3 Hi, I, + Hi, Hi, ,3) iC'^\S \ {b}) - C'^\S \ {c})) 



+ Hi, (5,3 , 5,, , i?Jt (5 \ {6}) - 5,3 , 5,, , i?f (5 \ {c}) ) 



4,4 



^3 h 
,4,4 



'3 h ^ 
4,4 



+ ,3 (5,, , 5,, , Bll (S \{b})- 61, , 5,, , (5 \ {c}) ) 



h h 
,4,4 



'2 U ^ 
4,4 



+ Hi, (5,, , 5,3 , 5^:1 (S \{b})- 61, , 5,3 , Bf'l (S \{c})) 



'h h 
,4,4 



^2/3^ 
4,4 



+ Hi, ,3 (5,, , 5,, , Bt;l (S \{b})- 61, , 5,, , 5^;!, (5 \{c})) 



+ Hi, (5,, , 5,3 , BfX {S \{b})- 61, , 5,3 , Bt;l {S \{c})) 



4,4 



+ Hi, (5,, , 5,, , BfX {S \{b})- 61, , 61, , (5 \{c})) 



4,4 



+ fe Si, h 5,3 b 5,4 b y4;^' ^ {S\{b})-5i,Ji,, S13 c Su c ^il I, ,3 ,4 



4,4 



(D.140) 



D.2 Relations between reduction coefficients 

The following relation is useful to cancel 1/B terms in the form factors for 5-point integrals: 

S;l {b^ - b,) - 6. - Svl) = (D.141) 

To prove eq. ()D.141|) . we introduce the auxiliary relation 



is^''^)u+s;i' + s;,'y, 



(D.142) 
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where y,,= ^ic 

ies\{c} 

To show relation ()D.142|) . let us assume the right-hand side is not zero, but some tensor ajic '■ 

- S^^^-^ji + S-^ji + S~\i yj, = aji, (D.143) 

Now, we compute 

''^^ <Ski ctjic = — Ski'S^'^^^^ji + Skj + Sckyjc 

les i€S\{c} 

-6j fc + fc = a k c 

-yjc + yjc = iik = c 
Since S is invertible, we must have ajic = 0. Summing eq. ()D.142|1 over / G S yields 

-bf + bj + bcyjc = (D.144) 

Now we multiply eq. ()D.142|) with be and eq. ()D.144|) with S~i^ and take the difference of the 

two resulting equations to obtain eq. ()D.141|) . 

Multiplying eq. ()D.141|) by A^*^ and summing over / G S leads to 

C^bf - b.C^^^l = C^e^b, - 6,C;, (D.145) 
Summing eq. ()D.145|1 over j in 5* \ {c} yields 

C^f, - 6, Vi^>: = C^f, i? - 6, (D. 146) 

and summing eq. ()D.145|1 . multiplied by A^^, over j in 5* \ {c}, we obtain 

v^'^K -\hc = - ^^^€1 ■ ^^•14'^) 

Further, one can show by direct calculation: 

5^6,(AL-"^D = l + Bml (D.148) 

i&S 

D.3 Special relations for = 5 

For N = 5, the external vectors form a basis of Minkowski space, such that the metric (in 4 
dimensions) can be expressed by the tensor H'^'^, constructed from external vectors only: 

(\ ^'^ 
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This fact implies the relation 



B 



ioT N = 5 



The proof is straightforward: 

Multiply the definition of H'^^ by A^^ AJ^^, we obtain 



E 



B 



jk 



^ja ■ \ l ^kb ■ \. 



Using now 
we get: 



a ' A j c ^ ( Sj c S j i -\- Sa i Sa i 



-A A — — A'^ 

ab '-^il '-^mn ~ o '-^mn y[4] '-^u ^mn 



2 — — iiLiL 2 

As the vectors A form a basis of Minkowski space for = 5, we conclude that 

1 
2 

On the other hand, using the definitions (see eqs. (j53|) and (j5T|l ^ 



''■ab — 9 y[4] 



^ ab 



g^-' + 2j2 S-\kKa^]b 
j,kes 



we see that T-f^ is equal to 



n 



keS 



^ 2 9[4] 



fJLV UU 



such that, together with eq. ()D.154j) . we must have 



D.4 Special relations for = 6 

B = j2h = 

i=l 

1=1 



^ for A^ = 5 . 



2 2\ 

mf 



i=l 



V A'' cS-^A"^. 



111/ 
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E Hexagon relations from helicity decomposition 



In the present appendix we provide very compact expressions for the reduction coefficients 
and Cf^ in the specific six-point case where the particles in the loop as well as all legs are 
massless. 

For massless gauge theory amplitudes usually the spinor helicity formalism |HH1IMIISH1IHB1 I57] 

is used for the treatment of vector bosons as well as for massless fermions. In this formalism 
the spinor degrees of freedom are systematically projected on helicity eigenstates. This is often 
essential if one wishes to write down the final result for an amplitude in a compact way. Now, 
one of the technical points which had emerged in the massless six-point calculations of jZU [Z2] 
is that the reduction coefficients bi can be written very compactly in terms of spinor traces, e.g. 



Here and in the following (12 • ■ ■) is a shorthand notation for tr {^^ ^2' ' ')■ The corresponding 
formulae for 62? • • • ? &6 are obtained by cyclic permutation. 

In the present appendix we first show that the helicity decomposition leads to further 
simplifications in eqs. ()E.161|1 . ()E.162|) and yields even more compact expressions. The two 
main identities responsible for the surprising simplifications in the Yukawa model calculation 
of jTTj become quite transparent in this approach. More surprising is perhaps that, even for 
the plain six-point determinant det(S, the introduction of chirality leads to a simple formula 
for its square root which cannot be written in terms of the Mandelstam variables alone. We 
present this further simplification of eqs. ()E.161|) . ()E.162|) in eqs. ()E.165|) . ()E.168|) below. The 
helicity representation also yields very compact expressions for the reduction coefficients C^^. 

We use here the helicity representation only as a means to get more compact formulae than 
those provided by the general method described in the main text for the fully massless = 6 
case. Before we proceed let us however stress that the helicity representation naturally leads 
to an alternative algorithm for the reduction of fully massless six-point rank r integrals, still 
into a linear combination of five-point rank r — 1 integrals, though decomposed on a different 
tensor basis. In this appendix we only sketch how this algorithm proceeds for the six point 
rank one tensor integrals. We then match the reduction thus obtained with the one derived in 
section which provides us with the above mentioned compact expressions for the reduction 
coefficients Cf„. 

Let us introduce the notation^ 



The product of two oppositely handed traces containing the same substring of some momenta 
Ph ■ ■ ■ Pii obeys the useful identity 

(+«i«2 ■ ■■iijij2 ■ ■ ■js){-iii2 ■ ■■iikik2 ■ ■ ■ h) 

^Our spinor helicity conventions follow |88| : in particular tr {-y^ijkl) — 4:ie{ijkl). 



detS 



(123456) (3456) - 2534^45556(6123) 
det5 

4S12S23S34S45S56S61 - (123456)^ . 



(E.161) 
(E.162) 
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= Si^i^Si^i^ ■ ■ ■ Si;_ii,((-)'jij2 ■ ■ -jskkikt^i ■ ■■kiii) (E.164) 
We can apply this to rewrite det S as follows: 

det5 = 4S12S23S34S45S56S61 - (123456)^ 

= 4(+123456)(-123456) - (^(+123456) + (-123456)) 

= -(^(+123456) - (-123456))^ = -(75123456)^ (E.165) 

Similarly, in eq. ()E.16H) the numerator of, say, 65 can be written as 

(123456) (1234) - 2si2S23S34(4561) = (^(+123456) + (-123456)) (^(+1234) + (-1234)) 

-2(+123456)(-1234) - 2(-123456)(+1234) 
= (^(+123456) - (-123456)) ((+1234) - (-1234)) 

(E.166) 

Combining this result with eq. ()E.165|) one factor of (+123456) — (—123456) cancels between 
numerator and denominator, and 65 becomes simply 

^ (+1234) - (-1234) _ (75I234) _ e(1234) 
^ (+123456) - (-123456) (75I23456) (75I23456) ^ ' ' 



By cyclicity. 



.((. + 2)(. + 3)(» + 4)(. + 5)) 
' (75 123456) ^ ^ 



The redundancy relation ^^^^ &i = (eq. ()D.157|) ) is straightforwardly checked by the chiral 
representation ()E.168|) for the reduction coefficients 6j, as it now reads: 

£(1234) - £(2345) + e(3456) - ^(4561) + ^(5612) - ^(6123) = (E.169) 

which is ensured trivially by momentum conservation^. 

Next, let us recall the two main identities used in [7T] for simplifying the calculation of the 
six-scalar Yukawa model amplitude: 

(3456) + (123456)6i = -2S34S45S5664 (E.170) 
(14) +62(1234) +65(4561) = (E.171) 

Eq. ()E.170|) is proved as follows: 

[(+3456) + (-3456)] [(+123456) - (-123456)] 



(3456) + (123456) 6 



^+123456) - (-123456) 



^It is not meant here that eq. (|D.157|I is the consequence of momentum conservation. On the contrary, it 
expresses the extra constraint that the Gram determinant vanishes for N = 6 |i46i i50j i.e. that any five of the 
six external momenta pi, i = 1, • • • , 6 are linearly dependent in a four- dimensional space-time. 
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[(+123456) + (-123456)] [(+3456) - (-3456)] 
(+123456) - (-123456) 
2(+345612)(-3456) - (-345612)(+3456) 
(+123456) - (-123456) 
(+1236) - (-1236) 
- 2S34S45S56 (+123456) - (-123456) 

= -2S34S45S56&4 (E.172) 

Here (IE.164jl was used in the third step. The proof of ()E.171|) is similar. 

The representation ()E.168|) for bi further suggests that, for the reduction of six-point tensor 
integrals, it might be natural to expand one of the numerator momenta in terms of vectors dual 
to the external momenta along the lines of [101. Let us define a dual basis {vi\v2\vl^\v^^}, 
i = 1, . . . , 6, for each of the six choices of four consecutive momenta^": 

= (-l)V(/i,z + 3,i + 4,z + 5), 

41 = (-l)V(^ + 2,/i,^ + 4,z + 5), 

41 = (-l)V(^ + 2,2 + 3,/i,2 + 5), 

41 = (-l)V(« + 2,z + 3,^ + 4,/i) 

For generic external momenta, each of these six bases can be used to expand the loop momentum 
k. Namely, further defining 



e'^^ = (-l)V(i + 2, i + 3, i + 4, i + 5) (E.173) 



one has 



Therefore, 



■ Pk = e^'^5k,,+a+i {k = t + 2,...,i + 5) (E.174) 



C = ^J24^'P^+^+^-^-- (1 = 1,. ..,6) (E.175) 

a=l 

The only caveat here is that we have used a projection of the momentum q^^ in the loop on 
four-dimensional space, while the loop integration requires dimensional continuation. However, 
for rank one integrals this is safe, since at one loop the index fi will be contracted with some 
external momentum or polarisation only. We will comment briefly on the case of higher ranks 
below. 

Using this equation together with eq. ()D.159|) and rewriting 

1 



Pi+a+l ' Qam 2 



Qi+a+1 ~ Qi+a ~ <Si+a+l,am + '^j+a.a™^ (E.176) 
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Note that contrary to j33| we use the pi rather than the n in building the dual basis vectors. 
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we obtain 

4 



i=l a=l 



Yl ^ Yl ^« ■ 

i=l a=l 

2^ 6 4 

(75I23456) ^ Yl \_^i+a+l ~ 9i+a ~ <Si+a+l,am + <Si+a,am "^a ''^ (E.177) 



Note that the e*^*-* from the numerator of bi in eq. ()E.168|) has cancelled against the same e'-'^ 
from the denominator of eg. ()E.175|1 : this is the main point of the algorithm. 
The Va'^ are not independent since 

= vf^ (E.178) 

and, by momentum conservation, also 

v^-vl^ = -{vf^-vt'''') (E.180) 

Remarkably, the three identities ()E.178|) . ()E.179|) . ()E.180|) together just imply that the 48 prod- 
ucts of qf's in flE.177|) cancel out in pairs. This leaves us with 

Using the decomposition ()E.181|) in the rank one massless six-point integral 1^'^ we obtain the 
following reduction to five-point rank zero integrals: 

We can write this result more compactly by introducing the convention that v^''^ = v^"^^ = 0, 
and defining 

d^:^" = v^-^-v^ (z,a = l,...,6) (E.183) 
In terms of the vectors da^^, eqs. ()E.178|) . ()E.179|) . ()E.180|) can be neatly combined to 



J 

We also note that 



df> = -4^ (z,j = l,...,6) (E.184) 



Y^a^' = = !,•••, 6) (E.185) 

a=l 

Y.r';d!f^ = ^^Sr4l (^ = l'---'6) (E.186) 
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Now, further defining 



'^l^ = Y.^^-^^'-a' (^,a^ = l,...,6) (E.187) 



a=l 

one has 



2i 

^6 ' ('^m! "S") 



rV^^r. ^{S\{i}) (E.188) 

(75I23456) 5 V \ 1 j; V ; 



Comparing eq. ()E.188|) with eq. of section for rank one - in which the first term on the 
r.h.s. vanishes as there is no higher dimensional integral for = 6 - we identify 

Cf„ = V''^ (E. 189) 

''^ (75 123456) ^ ' 

Let us add one comment on the application of this algorithm in the case of arbitrary rank. 
Applying the decomposition (jE. 181^1 to the six-point integral of arbitrary rank we obtain 



C'^--^^(ai,...,a.;^) = - (^^123456) ^ 4" '""-'"^(«2, ■ ■ ■ , a.; ^ \ {i}) (E.190) 



In the case of arbitrary rank the implied projection of the loop momentum on four- 
dimensional space requires a more careful consideration. This trick could fail if a contraction 
of two indices with a metric tensor g^^ occurs, as can be the case at ranks r > 2. However, 
if g^v is the full n-dimensional metric, then the integral breaks down to a lower rank integral 
anyway, while if g^^ is the (n — 4)-dimensional metric then g[n-4\ = g[n] — fi'[4] can be used before 
applying (IE.19()j) . 

In any case the index "/x" gets absorbed into the vectors (ii*^^, and ultimately will be 
contracted either (i) with an external momentum p^^, (ii) with an e - tensor, or (iii) with a 
polarisation vector ef^{pk,q)- Let us consider these cases in turn. 

(i) Contraction with external momenta produces e-tensors. These can, using momentum con- 
servation, be reduced to the six e^^^'s. 

(ii) Here we encounter the contraction of two e-tensors involving arbitrary momenta, 

(Pa -Pi Pa- Pj Pa " Pk 
Pb -pi Pb- Pj Pb ■ Pk 
Pc -Pi Pc- Pj Pc ■ Pk 

(abckji) — (abcijk) (E.191) 



(iii) As usual we take the reference momentum for si to be some other external momentum ps, 
so that 
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Then it is easily shown by using the Fierz identity that 



i \{ks) 



i+Hjk) 



i-Hjk) 



^l^,iJ>s)e{^i,i,i, k) 



2V2{ls)\-{kl) 
i \[ks]. 



i-Ujk) 




(E.193) 



2 V2[/s] L [kl] 



In writing eqs. ()E.193|) we have assumed i,k ^ I, which is not a restriction. 
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